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ABSTRACT 



Context. Mass-loss, occurring through radiation driven supersonic winds, is a key issue throughout the evolution of massive stars. 
Two outstanding problems are currently challenging the theory of radiation-driven winds: wind clumping and the weak-wind problem. 
Aims. We seek to obtain accurate mass-loss rates of OB stars at different evolutionary stages to constrain the impact of both problems 
in our current understanding of massive star winds. 

Methods. We perform a multi-wavelength quantitative analysis of a sample of ten Galactic OB-stars by means of the atmospheric 
code cmfgen, with special emphasis on the L-band window. A detailed investigation is carried out on the potential of Br ff and Pf r as 
mass-loss and clumping diagnostics. 

Results. For objects with dense winds, Br tt samples the intermediate wind while Pf y maps the inner one. In combination with other 
indicators (UV, H„, Br y ) these lines enable us to constrain the wind clumping structure and to obtain "true" mass-loss rates. For 
objects with weak winds, Br„ emerges as a reliable diagnostic tool to constrain M. The emission component at the line Doppler-core 
superimposed on the rather shallow Stark absorption wings reacts very sensitively to mass loss already at very low M values. On the 
other hand, the line wings display similar sensitivity to mass loss as H„, the classical optical mass loss diagnostics. 
Conclusions. Our investigation reveals the great diagnostic potential of L-band spectroscopy to derive clumping properties and mass- 
loss rates of hot star winds. We are confident that Br„ will become the primary diagnostic tool to measure very low mass-loss rates 
with unprecedented accuracy. 

Key words. Infrared: stars - stars: early-type - stars: winds, outflows - stars: mass-loss 
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1. Introduction 

In the last decade, massive stars (Mzams ~ 10M o ) have (re-) 
gained considerable interest among the astrophysical commu- 
nity, particularly because of their role in the development of the 
early Universe (e.g., its chemical evolution and re-ionization) 
and as (likely) progenitors of long gamma-ray bursters. Present 
effort concentrates on modeling various dynamical processes in 
the stellar interior and atmosphere (mass loss, rotation, magnetic 
fields, convection, and pulsation). Key in this regard is the mass 
loss that occurs through supersonic winds, which modifies evo- 
lutionary timescales, chemical profiles, surface abundances and 
luminosities. A well-known corollary in massive star physics is 
that a change of their mass-loss rates b y only a factor of two has 
a dramatic effect on their evolution (Mevnet et a 0l994l) . 

The winds from massive stars are described by the radiation- 
driven wind theory dCastor et al.l Il975l iFriend & Abbott! Il986l 
Pauldrach et al. 1986). Albeit its apparent success (e.g., 
Vink et al.ll2000l IPuls et al.ll2003l) , this theory is presently chal- 
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lenged by two outstanding problems (reviewed by iPuls et all 
2008), the clumping and the weak wind problem. 

The clumping problem. During recent years, various evidencd] 
has been accumulated that hot star wind are not smooth, but 
clumpy, i.e., that they consist of density inhomogeneities which 
redistribute the matter into clumps of enhanced density, embed- 
ded in an almost rarefied medium. 

Theoretically, the presence of such small-scale structurd3 has 
been expected since the first hydrodynamical wind simulations 
(lOwocki et al.lll988h . due to the presence of a strong instabil- 
ity inherent to radiative line-driving. This can lead to the devel- 
opment of strong reverse shocks, separating over-dense clumps 
from fast, low-density wind material. Interestingly, however, 
the column-depth averaged densities and velocities r emain very 



close to the predictions of stationary theory (see also Feldmeier 
1995. For more recent results, consult iRunacres & Owocki 



2002, 2005 (1-D) and lDessart & Owockill2003l 120051 (2-D)\M 
least for OB-stars, however, a direct, observational evidence in 
terms of line profile variability has been found only for two ob- 



1 For details, we refer to the proceedi ngs of the internatio nal work- 
shop on "Clumping in Hot Star Winds" (Ha mann et all 2008). 

2 not to be confused with large scale structure which is indicated by 
the ubiquitous presence of recurrent wind profile variability in the form 
of discrete absorption components (DACs, e .g.,|Prinia & Howarth [l986l . 
Kape r et alJl 19961, lLob el & Blomme 2008) and "modulation features" 
(e.g jFullerton et al.lll997fc 
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jects so far, the Of stars 4" Pup and HD 93 129A dEversberg et alJ 
119981 iLenine & Moffadl2008f) . 

Indirect evidence for small-scale clumping, on the other 
hand, is manifold, and is mostly based on the results from quan- 
titative spectroscopy, using NLTE model atmospheres. In order 
to treat wind-clumping in the present generation of atmospheric 
models, the standard assumption of the so-called "microclump- 
ing model" relates to the presence of optically thin clumps and 
a void inter-clump mediumQ A consistent treatment of the dis- 
turbed velocity field is still missing. The over-density (with re- 
spect to the average density) inside the clumps is described al- 
ternatively by a volume filling factor, / v , or a clumping factor, 
/ci > 1 , which in the case of a void interclump-medium, are 
related via / v = f~ x , The most important consequence of such 
a structure is that any mass-loss rate, M, derived from density- 
squared dependent diagnostics (H ff , Br ff or free-free radio emis- 
sion, involving recombination-based processes) using homoge- 
neous models needs to be s caled down by a factor of -Jfa. 

Ba sed on this approach. |Crowther et alJ d2002l) . lHillier et al] 
(l2Q03h and lBouret etal.1 (l2003l 120051) derived clumping factors 
of the order of 10. . . 50, with clumping starting at or close to 
the wind base. From these values, a reduction of (unclumped) 
mass-loss rates by f actors 3. . . 7 seems to be necessary (see also 
iRepolust et aTll2004l) . 

Even wor se, the analysis of the FUV Pv-lines by 



iFullerton et al. (2006) se ems to imply facto rs of 10 or even more 
(but see also Waldron & Cassinelh1 l2010l who argued that the 



ionization fractions of Pv could be seriously affected by X UV 
radiation). However, as suggested bv lOskinova et alJ (120071) . the 
analyses of such optically thick lines might require the consider- 
ation of wind "porosity", w hich reduces the effe ctive opacity at 
optically thick frequencies dOwocki et al.l l2.004). Moreover, the 
poros ity in velocity space (= "vorosity") might play a role as 
well (Owocki 2008). Consequently, the reduction of M as im- 
plied by th e work from Full e rton e t al. might be overestimated. 

Indeed, Sund qvist et alJ (1201 Oh . relaxing all the above stan- 
dard assumptions, showed that the microclumping approxima- 
tion is not a suitable assumption for UV resonance line for- 
mation under conditions prevailing in typical OB-star winds. 
These results a re supported for the case of B supergiants by 
Prinia & Massa] d2010h . who found that the observed profile- 
strength ratios of the individual components of UV resonance 
line doublets are inconsi stent with lines form ed in a "micro- 
clumped" wind (see also Sundqvist et al.l l201ll). Resonance to- 
gether with H„ line profiles as calculated by Sundqvist et al. 
d2010l 1201 ll) from 2/3D, stochastic wind models allowing for 
optically thick clumps (= "macroclumping") and vorosity ef- 
fects are compatible with mass-loss rates an order of magnitude 
higher than those derived from the same lines but using the mi- 
croclumping technique. 

Low mass-loss rates as implied by the latter models would 
have dramatic conseque nces for the evolution o f and feed-back 
from massive stars (cf. ISmith & Owockill2006l) . iHirschil d2008l) 
concluded that evolutionary models could "survive" with M re- 
ducti ons of at most a factor of ~2 in comparison to the rates 
from lVink et al.l ([2000) (which translate to an "allowed" reduc- 
tion of the empirical mass-loss rates of a factor of about four) 
whilst factors around 10 are strongly disfavored. Furthermore, 
such revisions would cast severe doubts on the theory of radia- 
tive driving, since the present agreement between observations 
and theory would break down completely. 



Hence, a reliable knowledge of the amount of clumping 
(quantified by the clumping-factor and its radial stratification) 
is crucial to constrain the "true" mass-loss rate of the star. Since, 
due to different oscillator strengths and cross-sections, the cor- 
responding formation depths vary from close to the base (Hq.) 
over intermediate regions (Br„, mid-IR continua) to the outer- 
most wind (radio), a consistent analysis of different diagnostic 
features will provide severe constraints on the run of f c \ and M 
itself. To this end, we have started a project to exploit these diag- 
nostics, by collecting the required data and analyzing them in a 
consistent way. First results with respect to constraints from the 
IR /mm/radio-contin uum combined with H a have been reported 
bv lPuls et al.l y006), in particular regarding the radial stratifica- 
tion of the clumping factor. They found that, at least in dense 
winds, clumping is stronger in the lower wind than in the outer 
part, by factors of 4. . . 6, and that unclumped mass-loss rates 
need to be reduced at least by factors 2. . . 3, in agreement with 
the results quoted above. 



The weak wind problem. From a detailed UV-analysis, 
iMartins et alJ (120041) showed that the mass-loss rates of young 
O-dwarfs (late spectral type) in N81 (SMC ) are significantly 
smaller than predicted theoretically (see also lBouret et alJ 2003 
for similar findings), even when relying on unclumped mod- 
els (the presence of clumping would increase the discrep- 
ancy). In the Galaxy, the same dilemma see ms to apply, 
articul arly for objects with log(L/L Q ) iS 5.2 dMartins et al.l 
2005b] iMarcolino etail 120091) . including the 09V standard 
10 Lac (|Herrero et al.ll2002l) . and maybe also r Sco (B0.2V, see 
IRepolust et al.l2005l) . pointing towards very low mass-loss rates, 
thus challenging our current understanding of radiation-driven 
winds. Note that most mass-loss rates for (other) dwarfs de- 
rived so far are only upper limits, due to the insensitivity of the 
usual mas s-loss estimator Hq, on (very) low mass-loss rates (see 
also Mokiem et al. I I2006I) . Present results based on UV studies 
may suffer from effects such as X-rays, advec tion or adiabatic 
cooling, as discussed by IMartins et alJ d2005bl see also iHillieri 
2008, lPuls et alj20"08l) . Consequently, a detailed investigation by 
means of sensitive mass-loss diagnostics in dwarfs over a larger 
sample is crucial to confirm their very weak-winded nature. 
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3 The importance of a low-density inter- clump medium for t he pro- 
duction of Ovi has been outlined already by Zsareo et al. (2008). 



In this paper, we intend to show that IR spectroscopy, in particu- 
lar in the L-band, is perfectly suited to investigate both problems, 
due to the extreme sensitivity of Br„ on mass-loss effects. 

i) For objects with large M, this line samples the intermedi- 
ate wind (because of the larger oscillator strength of Bi a 
compared to H ff and Br r ), thus enabling us to derive con- 
straints on the (local) clumping factor, and, in combination 
with other indicators (UV, Hq., Br r , Pf r ), to derive "true" 
mass-loss rates. We are aware that our UV-analysis might 
be hampered by macroclumping/vorosity effects (see above), 
but lacking suitable methods to include these effects into our 
NLTE treatment, we consider corresponding constraints as 
suggestive only. 

ii) For objects with very weak winds, Br„ provides not just up- 
per limits but reliable constraints on M. The relevance o f 
Br a has been pointed out already bv lAuer & Mihalasl (11969). 
who predicted that even for hydrostatic atmospheres the (nar- 
row) Doppler-cores should be in emission, superimposed on 
rather shallow Stark-wings. As we will show in the follow- 
ing, this emission component (and the line wings!) react sen- 
sitively on M, par ticularly for very weak winds (see also 
lNajarroetal]ll998l) . 
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Table 1. Sample stars and observing data in the L-band. "S" and "I" correspond to the IRTF/SpeX and VLT/ISAAC spectrograph, 
respectively, where HD 37128 and HD 37468 have been observed with both instruments. 





spectral 


sp. type 




obs. 


integr. 


previous investigations 


star 


type 


reference 


instr. 


date 


time 


Pv opt 


H/K IR/radio 


Cyg OB2 #7 


03 If 


MT91 


S 


10, 11 Sep 05 


360s, 400s 


3,6 


7 8 


CygOB2#8A 


05.5 1(f) 


MT91 


S 


11 Sep 05 


400s 


3,6 


7 8 


CygOB2#8C 


05 If 


MT91 


S 


10, 11 Sep 05 


360s, 400s 


- 3,6 


7 8 


HD 30614 (a Cam) 


09.5 la 


W72 


S 


11 Sep 06 


200s 


1 4,5 


7 8 


HD 36861 OlOri A) 


08 III((f)) 


W72 


I 


08 Jan 06 


816s(3.7yum ), 1224s(3.9jum ) 


1 4 


8 


HD 37128 (e Ori) 


BO la 


WF90 


I 


08 Jan 06 


102s(3.7/mi ), 204s(3.9^m ) 


2 


7 








S 


11 Sep 05 


80s 






HD 37468 (cr Ori) 


09.5 V 


CA71 


I 


08 Jan 06 


3068(3.7^ ), 612s(3.9/im ); 




7 








s 


11 Sep 05 


320s 






HD 66811(£ Pup) 


04 I(n)f 


W72 


I 


08 Jan 06 


204s(3.7//m ), 408s(3.9^m ) 


1 4,5 


7 8 


HD 76341 


091b 


M98 


I 


08 Jan 06 


826s(3.7^m ), 1656s(3.9//m ) 






HD 217086 


07 Vn 


W73 


s 


10, 11 Sep 05 


420s, 400s 


1 5,6 
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Spect ral references: C A71, IConti & Alschulej dl971l ) ; M98 , iMason et al.l dl998h : MT91, iMassev & Thompsor] d!99 lh : W72, IWalbornl d!972h ; 
W73. IWalbord dl973h : WF90. I Walbom & Fitzpatrickl dl99Ch . 

Previous inv estigations refer to ( 1) the analysis of the Pv/M 1118/ 28 doublet by Fullerton et al. (200(jj); (2) iKudritzki et al] d 19991 unblanketed 
analysis); (3) lHerrero et alJ l 2002): (4) the H g m ass-l oss analysis b y Markova et al. ( 2004) base d on s tellar parameters calibrated to t he results from 
optica l NLTE analyses bv lRepolust et alj d2004l) ; (5) lRepolust et alj d2004h : (6)lMokiem et alj d2005h : (7) the H/K-band analysis bv lRepolust et al.l 
(2005); and (8) the combined H^,, IR, mm and radio continuum analysis bv lPuls et al.l (12006). 



To accomplish our objective(s), we have performed a pi- 
lot study with the high resolution IR spectrograph ISAAC at- 
tached to the 8.1m Unit 1 telescope of the European Southern 
Observatory, Very Large Telescope (VLT), and the intermediate 
resolution spectrograph SpeX at the NASA Infrared Telescope 
Facility (IRTF), and secured high S/N L ( '-band spectra of se- 
lected Galactic OB stars. These spectra will be analyzed in 
the course of the present paper that is organized as follows: In 
Sect. |2] we describe our stellar sample, the observations and the 
data reduction. Sect. [3] summarizes the relevant features of the 
used NLTE model atmosphere, and describes our implementa- 
tion of clumping. In Sect.|4]we concentrate on those objects with 
dense winds within our sample, and investigate their clumping 
properties by combining the L-band spectra with other diagnos- 
tics. The complementary objects with thin winds are considered 
in Sect. [5] after the principal line formation mechanism of Br„in 
thin winds has been discussed, and additional problems have 
been illuminated. In Sect. [6j finally, we discuss our results (par- 
ticularly with respect to wind-momentum rates), and summarize 
our results and present our future perspectives in Sect. [7] 



2. The stellar sample, observations and data 
reduction 

The specific targets comprise a subsample of the northern 
(SpeX) and southern (ISAAC) objects of the large sample of 
Galactic OB-stars which has been observe d and analyzed in the 
optical (see Table[TJl and in the H/K-band (Repo iust et all 2005. 
based on the material presented bv lHanson et alJ 120051) . Our sub- 
sample covers mostly supergiants from 03 to B0, and has been 
augmented by one weak wind candidate (HD 37468, 09.5V) 
and three comparison objects (HD 217086, 07Vn; HD 36861, 
08III((f)), HD 76341, 091b) which should behave as theoreti- 
cally predicted. Note that for all objects UV archival data are 
available as well, that four of them have been analyzed with 
respect to Pv (cf. Sect. [TJ), and that all O-supergiants plus the 
08 giant have been investigated in H a and the IR, mm- and 



radio-continuum by iPuls et all d200 6) regarding their clumping 
properties. Two of our objects have been observed by SpeX and 
ISAAC, to enable a comparison of the data obtained by both 
instruments. Table [TJ summarizes our target list, important ob- 
serving data and the previous investigations of our objects in the 
individual wavelength bands. 

2.1. The IRTF/SPEX Sample of Stars 

Our first spectra were obtained in Septemb er, 2005 at the 3.0 m 
IRTF on Mauna Kea, Hawaii, using SpeX (iRavneret al.l l2003). 
Spex is a medium-resolution, 0.8 to 5.4 micron spectrograph 
built at the Institute for Astronomy (IfA) and is available for 
general use to the public through qualified time on the IRTF. 
Two nights were granted to our program. While the second night 
was good, the first night experienced intermittent cirrus, some- 
times seriously reducing the signal in the instrument. Our spectra 
were taken in the cross-dispersed mode of SpeX, providing full 
spectral coverage from 0.8 to 5.5//m with a resolution of A/6A ~ 
2000 with the narrowest slit (0.3")- Luckily, the seeing was never 
greater than about 0.8" and was typically closer to 0.6". About 
every hour during the night, wavelength arcs and flat fields were 
taken. Spectra were also obtained of telluric standards. These 
telluric standards were selected to share the same airmass and 
approximate sky location as the target stars, are observed to con- 
siderably higher counts, and were chosen to be early A-dwarfs 
with low vsim. 

A clear advantage to using the IRTF/SpeX system in cross- 
dispersed mode is the ease with which the spectra can be 
reduced. A sophisticated IDL-ba sed reduction pack age called 
SpeXtool has been developed by Cushi ng et al.l ([2004) which in- 
corporates calibration frames and readily produces final wave- 
length and flux calibrated spectra. While the data is cross- 
dispersed, there is sufficient room in the 15" slit to allow two 
uniquely-observed positions. This provides for a traditional A 
position - B position for background subtraction. SpeXtool also 
takes advantage of the telluric correction methods recently de- 
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veloped bv lVacca et al.l (l2003h which includes a high-resolution 
model of Vega. Because Vega has an extremely low v sin i, the 
A-dwarfs selected for tellurics in our study (HD 219190 and 
HD 33654) were also chosen to have very low v sin to ensure 
the be st match in profile shape. The Appendix of lHanson et alj 
(2005) provides graphic evidence of why the vsin; match be- 
tween model and star is so critical for very high signal-to-noise 
spectroscopic work like this. 

We were particularly lucky for this program as Vega was ob- 
servable at the start of the night during our run. We took full 
advantage of this for the purpose of carefully confirming the in- 
tegrity of the high signal-to-noise telluric spectra derived from 
all the A-dwarf standard stars observed throughout the night 
compared against the theoretically-perfect telluric spectrum de- 
termined from our direct Vega observations and derived by the 
Vega models provided in the SpeXtool package. 

While spectra were obtained throughout the full spectral 
range of 0.8 to 5.4 microns with SpeX, we are presenting just 
the two (most interesting) narrow spectral regions centered at 
Pf r and Bi a in Figure Q] 

2.2. The VLT/ISAAC Sample of Stars 

We were granted one night, in visit or mode, on 8 Januar y 
2006 on VLT1 (Antu) with ISAAC (Mo orwood etail Il998l) . 
The weather was at times marginal, with highly variable see- 
ing and sometimes cloud cover too dense to observe. To achieve 
the highest resolution, again we had the slit set to 0.3". This 
proved extremely challenging when the seeing dropped below 
2.0". However, during that single night, we also experienced a 
few extended moments of reasonable weather and seeing. By 
sticking to the brightest sources in our sample we were able to 
observe several stars with sufficiently high count rates, even at 
4.05/im, to achieve the high signal-to-noise needed for our anal- 
ysis. 

For reduction, we closely adhered to the advise fou nd in 
the ISAAC Data Reduction Guide 1.5 (lAmico et al.ll2002l) and 
updates found on the ESO ISAAC website. An all encompass- 
ing reduction package is not available for the ISAAC instrument 
such as is available with SpeX. We used fits manipulation rou- 
tines available from the IRAF0 software package. The reduction 
starts with a simple ESO provided configuration file to remove 
electrical ghosts (provided in the 'eclipse' package). From there 
reduction steps involved dark subtraction, linearity corrections 
and flat fielding, all accomplished using scripts written in IRAF. 
ISAAC has ample slit length for multiple positions in the slit. 
However, we observed all of our targets and our telluric stan- 
dards in the near exact same two positions in the slit, noted as 
position A and position B. 

The ISAAC instrument shows a pronounced curvature and 
distortion with wavelength on the array. Perhaps more impor- 
tantly, the shape of this curvature is also a function of position 
along the slit. Before extracting a 1-D spectrum from the 2-D 
image, this needs to be corrected. ESO has provided scripts in 
the eclipse package that use the arc images taken throughout the 
night to create a distortion correction that can be applied to the 
2-D images. For typical applications, this allows stars observed 
anywhere in the slit, to line up properly in wavelength space. 



4 IRAF is distributed by the National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities for Research 
in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 



However, these corrections will not be good enough for our ob- 
servations. 

This is the reason only two positions were used in the slit. 
Two positions are needed to remove the background in the 2-D 
images. The distortion in the 2-D images was corrected using 
the ESO eclipse program for this purpose. But we were care- 
ful to use a single solution applied to all 2-D images taken dur- 
ing the night, so there was no introduction of very small, but 
differing wavelength solutions. Once the 2-D image was fully 
processed and had been converted to a 1-D spectrum, the two 
slit positions, A and B, were never used together in processing 
again. Target star A slit positions were only processed with tel- 
luric star A slit position and the two slit positions, A and B, were 
kept separate in processing in much the same way two differ- 
ing grating tilts would proceed separately. At our resolution and 
signal-to-noise, subtle misalignment of grating solutions in this 
wavelength regime, because of the numerous, very deep and nar- 
row telluric features would create strong beat phenomena when 
1-D spectra are divided from each other in later processing. 

To remove telluric absorption, our strategy is one of 
bootstrapping all telluric observations off each other (see 
lHanson et al.ll2005l) to come up with a consistent set of telluric- 
free spectra. We started with a synthetic spectrum of the Earth's 
atmosphe re for the airmass covering our observations using 
ATRAN (lLordlll992l) . This model telluric spectrum was used 
to divide out (remove) the telluric component in all of our tel- 
luric A-stars to first order. We then fit the remaining hydrogen 
lines in the A-stars. We also went through this exercise using 
several of the OB-stars. While we did not derive the hydrogen 
profiles of these OB stars in this manner, the relatively narrow 
width of the OB hydrogen lines allowed us to use their spectra 
to constrain the very broad wing component of the A-star and to 
ensure a proper continuum for the hydrogen line fit in the A-star. 
Then we returned to the original full A-star spectrum, removed 
the fitted hydrogen lines for that star and created what was the 
best estimate for the telluric features through that spectral range 
towards that star. In this way, the telluric spectrum was individ- 
ually solved for numerous A-star sight-lines with a similar air- 
mass range. These can be checked against each other and then 
averaged to reduce any possible errors introduced in the hydro- 
gen line fits from a single A-star. From this, a final few telluric 
spectra as a function of airmass during the night was derived, 
and the appropriate telluric spectrum could be removed from the 
OB target star. A rough, but independent check was made by di- 
viding our raw OB star spectra by an ATRAN derived telluric 
spectrum to ensure there were no obvious mistakes introduced 
in our method above. 

The validity of the reduction procedure for each dataset is 
confirmed in Fig. [2] where a comparison between IRTF/SpeX 
and VLT/ISAAC L-band spectra for the two objects (HD 37468 
and e Ori) observed in both runs is presented. In Fig. [2] the 
VLT/ISAAC spectra have been reduced to the resolution of the 
IRTF/SpeX instrument (note the larger S/N ratio of the former). 
The excellent agreement between both datasets clearly supports 
the different data reduction procedures employed for each run. 
Given the higher resolution of the VLT/ISAAC data, we made 
use of these observations for our quantitative analysis. 

3. IR diagnostics 

To model the infrared spectra of our sample of objects we 
have utilized cmfgen, which is an iterative, non-LTE, line- 
bl anketed model atmosphe re/spectrum synthesis code developed 
bv lHiliier& Milled (119981) . It solves the radiative transfer equa- 
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CygOB2-8A 



HD217086 




HD37468 




eps Ori 



3.72 3.74 3.76 
wavelength [micron] 




0.2 



03If* (IRTF/SPEX 




04I(n)f (VLT/ISAAC) 



05.51(f) (IRTF/SPEX) 



07Vn (IRTF/SPEX) 



08111(f) (VLT/ISAAC) 




09.51a (IRTF/SPEX) 



09.5V (IRTF/SPEX) 



09.5V (VLT/ISAAC) 




BOIa (VLT/ISAAC) 



4.00 4.02 4.04 4.06 4.08 
wavelength [micron] 



Fig. 1. VLT/ISAAC and IRTF/SpeX L-band spectra for our sample stars in the two important regions centered around Pf y (left) and 
Br ff (right). Note the Hei3.703yum line which is present at later spectral types. Two objects (HD 37468 and e Ori) have been observed 
by both instruments. 
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Fig. 2. Comparison between IRTF/SpeX (solid) and 
VLT/ISAAC (dashed) L-band spectra for the two objects 
(HD 37468 and e Ori) observed in both runs. The VLT/ISAAC 
spectra have been reduced to the resolution of the IRTF/SpeX 
instrument. Note the excellent agreement between both datasets. 

tion in the co-moving frame and in spherical geometry for the 
expanding atmospheres of early-type stars. The model is pre- 
scribed by the stellar radius, the stellar luminosity, L„, the 
mass-loss rate, M, the velocity field, v(r) (defined by v m and /3), 
the volume filling factor characterizing the clum ping of the stel- 
lar win d, and elemental abundances. Following Hillier & Miller 
(1998, see also iPauldrach et al.ll 19941) we include X-rays char- 
acterizing the X-ray emissivity in the wind by two different 
shock temperatures, velocitie s and filling factors. We refer to 
iHillier & Miiierl dl998l 1 19991) for a detailed discussion of the 
code. 

Given the large range of stellar parameters and the variety 
of luminosity classes covered by our O-star sample, the location 
of the tr oss = 2/3 radius for these objects will vary from being 
placed in the deep hydrostatic layers (dwarfs) up to the upper 
layers where the wind takes off (supergiants). Noting that not 
only the IR lines but also the IR continuum, through bound-free 
and free-free processes, will have different formation depths as a 
function of wavelength, the role of the hydrostatic structure and 
the transition region between photosphere and wind becomes 
crucial to interpret the stellar spectra. 

With this in mind, we computed CMFGEN models with a 
photospheric structu re modified following the approach from 
ISantolava-Rev et ail dl997 ) 0, smoothly connected to a beta ve- 
locity law. In our approach the Rosseland mean from the orig- 
inal formulation was replaced by the more appropriate flux- 
weighted mean. Several comparisons using " exact" photospheric 
structures from tlusty dHubenv & Lanzl 19951) showed excellent 
agreement with our method. Likewise, model atoms were ex- 
panded and optimized to make use of the IR metal lines arising 
from high lying levels. 

5 In this approach we compute the density from the hydrostatic equa- 
tion and the velocity from the continuity equation. 



To investigate the role of clumping, we follow the conven- 
tional approach of microclumping and assuming a void inter- 
clump matter. In this case, and as already outlined in Sect. [1] 
the volume filling factor is just the inverse of the clumping 
factor, / v = / d ', and the clumping factor itself quantifies the 
overdensity of the clumps with respect to the averaged density, 
(p) = M /(4nr 2 v(r)). Moreover, the radial stratification of the 
volume filling factor i s desc ribed by the clumping law intro- 
duced by iNajOToetin (Eool): 

_ rlf) _ vo„-i<r) 

Mr) = CLi + (1 - CLi)e + (CL 4 - CL x )e =5 (1) 

where CLi and CL4 are volume filling factors and CL2 and CL3 
are velocity terms defining locations in the stellar wind where 
the clumping structure changes. CLi sets the maximum degree 
of clumping reached in the stellar wind (provided CL4 > CLi) 
while CL2 determines the velocity of the onset of clumping. CL3 
and CL4 control the clumping structure in the outer wind. From 
EqQ]we note that as the wind velocity approaches v«,, so that 
(vco - v(r)) < CL3, clumping starts to migrate from CLi to- 
wards CL4. If CL4 is set to unity, the wind will be unclumped in 
the outermost regi on. Such behavior wa s already suggest ed by 
Nugi s et al.l (119 98) and was utilized by iFiger et al.l (120021) and 
iNajarro et al.l (120041) for t he analysi s of th e WNL stars in the 
Arches Cluster. Recently, iPuls et af] J2006) have found a sim- 
ilar behavior from H a and IR/mm/radio studies for OB stars 
with dense winds. Furthermore, our clumping parametrization 
seems t o follow well the resu l ts fr om hydrodynamical calcula- 
tions by iRunacres" & Owockil (120021) . From Eq. Q] we note that 
if the term including CL3 and CL4 is neglected or if C L3 — > 0, 
we recover the law proposed by Hillier & Milled d 1999b . For the 
present study (except for ( Pup, see below), we have set CL4=1, 
i.e., the outer wind regions are assumed to be unclumped. 

Observational constraints are set by the L-band spectra pre- 
sented above and UV, high-resolution optical and H and A'-band 
spectra collected by our group as well as by optical, IR and radio 
continuum measurements from literature/archival data. The indi- 
vidual sources are quoted in the corresponding figure captions. 
In this paper we concentrate on the strong diagnostic potential 
provided by the infrared K- and L-bands to determine mass-loss 
rates and trace wind clumping as an alternative to H ff , the classi- 
cal mass loss indicator. Thus, we defer a detailed full wavelength 
analysis and discussion of these objects to a forthcoming paper. 

Table|2]displays the stellar parameters obtained for our sam- 
ple, whereas a detailed comparison with results from previous 
investigations is provided in Appendix lAl We obtain uncertain- 
ties of ~1000 K for the effective temperature of the objects (see 
Appendix [A] for a thorough discussion) while typical errors of 
0.1 dex are estimated for logg. For objects with dense wind, we 
estimate our mass-loss accuracy to be better than 25%, with the 
corresponding M//" 5 = const scaling for the error on the clump- 
ing factor. For objects with thin winds we consider 0.5 dex as a 
conservative error on the mass-loss rate estimate (see Sect. 16. lb . 

Projected rotational speeds, v sin i, have b een derived via the 
Fourie r-transform technique as developed by Simon-Diaz et al.l 
(2006) (based on the original method proposed by lGravill973l 
1975), applied to weak metal lines (and partly Hei-lines) avail- 
able in the spectra. The remaining discrepancies between syn- 
thetic and observed lineal were "cured" by convolving the spec- 
tra with an additional, radial-tangential macro-turbulent velocity 
distribution (entry v macro in Table [2J. The derived values are of 

6 Such di screpa n cies were alrea d y dete ct ed and descr i bed b y 
i Rosenhald d 19701). IConti & Ebbetsl dl977h . iLennon et all j!993l) . 
lHowarthetalJjl997h . 
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Table 2. Stellar and wind parameters as adopted/derived in the present analysis. T e ff in kK, R* in R , all velocities in km s" 1 , M 
in Kr 6 M /yr and modified wind-momentum rate, D m om = MvodR+fRof 5 , in cgs. The gravitational acceleration, \ogg, is the 
effective one, i.e., not corrected for centrifugal forces. The volume filling factor, / v , corresponds to the parameter CLi in Eq.[TJand 
describes the maximum degree of clumping reached in the stellar wind. However, since our analysis bases on stratified clumping 
factors, the usual scaling of M oc Jfc does not or only approximately apply in most cases. CL2 (km s _1 ) indicates the onset 
of clumping in the wind. The horizontal line separates objects displaying significant wind emission in H„, i.e., dense winds (see 
Sect.|4j>, from objects with a pure absorption H„ profile, i.e., thin winds (Sect. [5J. 
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similar order as found in alterna t ive investigations ([Rvans et ail 
20021 ISimon-DiazetalJ 120061 Isimon-Draz & Herrerd 120071 
Lefever et alj 120071 iMarkova & Pulsl 120081) . indicating highly 



supersonic speeds in photospheric regions tha t would be diffi- 
cult to exp lain. Recently, however, lAerts et all (I2009L see also 
iLucvl IT976h interpreted such extra-broadening in terms of col- 
lective effects from hundreds of non-radial gravity-mode oscilla- 
tions, where the individual amplitudes remain sub-sonic. QThey 
pointed out that the rotational velocity could be seriously under- 
estimated whenever the line profiles are fitted assuming a macro- 
turbulent velocity rather than an appropriate expression for the 
pulsational velocities, or if a Fourier technique is applied to in- 
fer the rotational velocity. If this were true, our values for v sin i 
(and also those from the quoted investigations) would provide 
only lower limits. For the present investigation, however, this is 
of minor concern, since our main interest is to obtain a correct 
shape of the profiles (irrespective of the responsible process), to 
enable meaningful fits. 

4. Objects with dense winds, constraints on the 
clumping factor 

In this section we discuss our results for the objects of our 
sample displaying dense winds, an d compare them wit h pre- 
vious studies carried out at optica l (|Repolust et al. 2004, here- 
after REP04: IMokiem et alJ I2005L MOK05) and near-infrared 
wavelengths (Repo lust et al.ll2005l REP05) as well a s the com- 
bined H„/IR/mm/radio analysis by iPuls et ail (120061 hereafter 
PUL06). For further details, see Appendix A. 
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Fig. 3. Model fits to IR and H ff observations of CygOB2 #7 
(solid, black). Two different models corresponding to clumping 
values of / v =0.03 (dashed, red) and / v =0.4 (dashed-dotted, blue) 
are displayed (see text). The H a spectrum was obtained with 
ISIS at the William Herschel Telescope on La Palma (kindly pro- 
vided by A. Herrero). 



Cyg OB2 #7. Our derived main stellar parameters for 
Cyg OB2 #7 (s ee Table |2| agree very well with those obtained 
from optical dMOK05l) and // and ST-band (REP05) analy- 
ses with respect to the effective temperatures (within less than 
1000 K). This is a very encouraging result since our temperature 
determination rel ies not only on the Hen/Hei equilibrium (e.g., 
lMOK05LlREP05h but also on the Nv/Niv/Nm equilibria. Our re- 
sults confirm the consistency between both criteria and the valid- 

7 First ob servational evidence in support of this scenario has been 
provided by Simon- DfazetalJd2010h . 



ity of the weak Hei lines as diagnostics in this high temperature 
regime. Likewise, our logg va lue lies in between the ones de- 
rived by MOKf3and|REP05, while our unclumped mass-loss 
rates (~ 7.8 -10" 6 M o /yr) are roughly 30% lower. We attribute 
this disc repancy to our lower He abundance (0.13 vs 0.21 in 
IMOK05I) and higher /3 (1.0 vs 0.8). Interestingly, our models fa- 
vor a large clumping starting relatively close to the photosphere 
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Fig. 4. Model fits (dashed) to IR and optical observations of ( Pup (solid). Optical data were retrieved from the ESO archive, 
programme 266.D-5655(A) 



to provide consistent simultaneous fits to the UV0, optical and IR 
observations of this object (see Fig. [3}. However, such a strong 
clumping - via the corresponding lower mass-loss rate - tends to 
produce much too deep cores in the optical H and Hen lines. For 
comparison, a less clumped model (/ v =0.4), better matching the 
optical lines, is also displayed in Fig. [3] Note, however, that such 
a model leads to severe mismatches with the L-band and UV (not 
displayed here) spectra. From Fig. [3] we see that in the case of 
strong winds Br y and Pf r provide stronger response to clump- 
ing than Br^. lPULOa also found a strong clumping in the inner 
wind of this object. Their average clumping factors lie between 
the ones we obtain in our IR and optical analysis. We would like 
to stress that while the optical and IR spectra of Cyg OB2 #7 
provide strong constraints on CLj and CL2, the UV spectra and 
submillimeter and radio observations constitute crucial diagnos- 
tics to det ermine CL4 and CL 3 . Indeed, our UV and submillime- 
ter data (Naj arro et al.l l2008) support the presence of constant 
clumping, at least up to mid-outer wind regions where the mil- 
limeter continua of Cyg OB2 #7 are formed. However, radio ob- 
servations bv lPULQq show that clumping may start to vanish at 
the outermost regions of the stellar wind (note that radio con- 
tinua form at much larger radii). The expected emission of our 
models with constant clumping severel y overes timate the upper 
limits provided by the observations by iPULOol of Cyg OB2 #7. 



8 remember our caveat regarding the impact of macroclumping on 
UV resonance lines, as stated in Sect.fT] 



This demonstrates the need of multiwavelength observations to 
constrain the run of the clumping structure. 



£ Puppis. Our derived main ste llar parame ters agree fairly well 
with those presented by REP04 and REP05. We find, however, a 
slightly lower He abundance and a much higher (by 50%) "un- 
clumped" mass-loss rate. We attribute this discrepancy in M to 
our clumping parametrization for this object. Given the large 
number of spectroscopic and continuum constraints at nearly 
all wavelengths available for £ Pup, we performed a detailed 
clumping study aiming to constrain as accurately as possible 
the run of the clumping factor th roughou t the wind and com- 
pare it with recent results from |PUL06l Thus, we made use 
of all the clumping parameters presented in EqQ] and obtained 
CLj=0.03, CL 2 =180., CL 3 =600 and CL 4 =0.17 as a best fit 
(see Fig. [T8l ). With this parametrization, maximum clumping 
(/v=0.03) is reached only in a very narrow, AR ^0.1 R*, re- 
gion around r=1.5 R t . Thus, the classical M/f^' 5 = const scaling 
for constant clumping does not hold and causes the above dis- 
crepancy regarding unclumped M values. Our best model (see 
Fig|4]i is able to reproduce satisfactorily not only the L- and K- 
band spectra but also the optical lines and the full UV to radio 
energy distribution of the object. We stress the almost perfect fit 
reached for H ff . That quality of fit can only be achieved for mod- 
els with stratified clumping as otherwise the observed absorption 
and emissio n components cannot be fitted simultaneously (see 
also lPUL06h . We note, however, that since we aimed at a com- 
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tors (~5 for Hem recombining at v=0.86vooand ~11 for a wind 
with completely ionized He, when assuming an unclumped outer 
wind) are reached in the same wind region (i.e. around r= 1.5/?*) 
as in our models. When scaled to a similar outer clumping fac- 
tor as derived in this work, the agreement is even more striking 
(see Fig.fTSl. Thus, both studies reach the same conclusions con- 
cerning the run of the clumping factor. As for Cyg OB2 #7, our 
investigation shows the high sensitivity of Br r and Pf r to clump- 
ing (see Fig. |4]i and the enormous potential of IR spectroscopy 
to constrain the structure of stellar winds. 



Cyg OB2 #8C. The derived stellar p arameters differ consider- 
ably from those obtained by MOK05. Our temperature is more 
than 4000 K lower while our gravity lies only 0.12 dex be- 
low. We are confident in our temperature determination since, 
as discussed above, we make use of both helium and nitrogen 
ionization equilibria (see Fig|5]l. In fact, our model reproduces 
satisfactorily the optical and IR spectra of this object. We sug- 
gest that th e discrep ancy is closely related to the lower vsin/ 
derived by 1MOK05I (v sin/=145km s vs. our 175km s + 
90 km s^'macroturbulence) and the strong reaction of Hei 4471 
in this parameter domain. The difference in derived v sin i values 
can be attributed to the fact that our spectra are of higher quality 
than those used by MOK05. We stress that our best compromise 
solution underestimates the emission core of Pf r . 

Compared to the other strong wind objects discussed in this 
section, this object requires a lesser degree of clumping. Further, 
our models imply an onset of clumping which is located at larger 
velocities than for the rest of our supergiants. IPUL06I also de- 
rived a low degree of clumping. However, a detailed comparison 
with their results is not possible as their analysis of this object 
remained rather unconstrained due to the lack of sufficient flux 
measurements. 



Cyg OB2 #8A. Unlike for the case of Cyg OB2 #8C, our anal- 
ysis of Cyg OB2 #8 A yi elds exce l lent agreement with the stellar 
parameters obtained by MOK05. de Becke r et alj (120041) report 
this object to be a 061 + 05.5III binary system. This can be 
clearly inferred from the absolute magnitude of the system dis- 
played in Tabl^l] In fact, no single best fit could be obtained to fit 
simultaneously the optical and IR observations of Cyg OB2 #8A 
(see Fig|6]l. Our preferred model, which reproduces better the IR 
spectra, is characterized by strong clumping (/ v =0.01) that leads, 
as in the aforementioned case of Cyg OB2 #7, to somewhat too 
strong absorption cores in the optical hydrogen lines. We note 
that the inclusion of clumping nicely removes the discrepancy in 
Br r found by REP05. On the other hand, a model with / v =0.1 
tuned to optimize the optical (see Figj6j severely underesti mates 
the emission components in the IR lines. The analysis of |PUL06l 
was hindered by the non-thermal nature of CygOB2 #8A's ra- 
dio emission. Nevertheless, their combined H„ + IR photome- 
try analysis yields a wind structure which is significantly less 
clumped than inferred in this paper. Both results could be only 
reconciled if the outer wind would be strongly clumped. 



promise solution, our model is unable to fit the blue wing of Br r . 
Thus, while a different clumping would yield a better fit to Br r , 
it would also worsen significantly the model fits to other lines. 

In their analysis of the clumping structure of O-star winds, 
IPUL06I pointed out that the derived clumping values for f Pup 
are strongly dependent on the assumptions regarding the run of 
the He ionization. Interestingly, their maximum clumping fac- 



a Cam. We fi nd excel lent agreement with the stellar parame- 
ters derived bv lREP05L A very high degree of clumping, start- 
ing very close to the photosphere, is required to match Pf y . In 
fact, this line turns into the most sensitive clumping diagnostic 
at the base of t he wind in O supergiants. Our result agrees qual- 
itatively with lPUL06f s conclusions on the run of the clumping 
factor. They found a moderate degree of clumping in the inner 
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Fig. 7. Model fits (da shed) to IR and op tical observations of a Cam (solid). Optical data are from the Indo-U.S. Library of Coude 
Feed Stellar Spectra dValdes et alj |2004). 



and mid wind regions. Figure[7]shows that our model can repro- 
duce satisfactorily the IR and optical spectra of HD 30614. 

e Qri. Several spe ctroscopic studi es of e O ri using cmfgen (e.g. 
ISearle et alJ l2008i and fastwind (REP05), yielding similar pa- 
rameters, have recently appeared in the literature. While the cm- 
fgen study made use of UV and optical data, the fastwind one 
used in frared spectroscopic observations alone. Interestingly, 
REP05 could derive only an upper limit on the effective tem- 
perature due to the absence of Hen IR diagnostic lines. Our re- 
sults, arising from a full UV to IR investigation and displayed 
in Tabl^2] revise down the effective temperature by roughly 
1000 K. A striking result revealed by Fig{8] is the requirement 
of strong clumping to match the IR spectra. Figure [8] demon- 
strates the failure of an unclumped (/ v = 1) wind to reproduce 
the L-band lines. As for most of the previous stars, such strong 
clumping leads to overestimated HI and Hei line cores. In the 
case of e Ori, however, this mismatch could be also due to the 
intrinsic line profile variability, as the IR and optical observa- 
tions where taken at different epochs. 



5. Objects with thin winds 

HD 21 7086. Our derived stellar parameters for this fast ro tating 
dw arf are in excellent agreement with those obtained by REP04 
an d lREP05l bv means of optical and IR spectroscopy. While their 
IR analysis could provide an upper limit to the mass-loss rate be- 



ing a factor of two lower than the optical one, no firm determi- 
nation of this parameter could be assessed. Our L-band data (see 
Fig. [9j clearly show how Br ff constitutes a much more powerful 
M diagnostic for low density winds than the previously used H a 
and Br r lines. Our UV through IR study yields a less clumped 
wind than for the case of the supergiants. Without the UV we 
could not have broken the M- clumping degeneracy though. 

HD 36861. This object displays the strongest wind within our 
sample of stars with low density winds, as to be expected from 
its 08 III((f)) spectral type classification. Ou r derived effective 
temperature is hotter tha n the one adopted by|PU L06 (based on 
the calibrations used by Markova et al. 2004J). Interestingly, no 
clumping is required by our models to match the IR and optical 
lines. This result is cons istent with one of the two possible so- 
lutions found by PUL06. Figure. [TOl (upper two panels) shows 
the excellent agreement of our models with the observations. 
Exceptions are two of t he optical Hei single ts (related to the so- 
called singlet problem. iNaiarro et al.l l2006) and the He compo- 
nents of the IR H lines. Our models (see Fig. [10] upper left) tend 
to show the Hei components in emission].] Since this discrep- 
ancy appears only in objects with low density winds where the 
Hei hydrogenic components form at higher densities, we suggest 



9 As a guideline, we plot as well the line profiles for Br„ and Pf r 
computed without Hei components (dashed-dotted) and with no He at 
all (long dashed). 
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Fig. 9. Model fits (dashed) to IR and H ff observations of 
HD 217086 (solid). H a observations as in Fig|3] 



that realistic broadening functions should be developed and used 
to replace the assumed pure Doppler profiles. 

HD 76341. As for HD 36861, no clumping is required to re- 
produce optical and IR (also UV, not shown here) spectra of 
HD 76341 (Figure [10] lower two panels). Again we note the 



Fig. 11. Model fits (dashed) to IR and H a observations of 
HD 37468 (solid). The H„ profile was obtained with IDS at the 
INTTelescope onLa Palma (kindly provided by S. Simon-Diaz). 



problem with the Hei components in Br„, but also the enormous 
potential of this line to determine mass-loss rate in thin winds 
when H ff struggles to react to changes in M. 

HD 37468. This is the object with the thinnest wind in our sam- 
ple. Our analysis has made use of available UV and optical spec- 
tra as well. Once more, no clumping is required and we obtain 
M- 2 •lO _lo M /yr as our current best estimate (see Fig. fTTT i. 
The need of correct broadening functions for the Hei lines is 
again evident in the Br ff complex. Even though our M determi- 
nation appears perfect, we stress that in this regime of extremely 
low mass-loss rates the resulting synthetic Br„. profile can be 
very sensitive to the data set used for the hydrogen collisional 
bound-bound processes (see below). 

In the following, we will discuss the formation of the specific 
shape of the Br„ profile in these very thin winds in considerable 
detail. 



5.1. Theoretical considerations 



As has been extensive l y discu ssed by Mihal a s| (ll978l) . Kudritzkil 
(I1979I). iNaiarroetal l (1 1998b . iPrzvbilla & Butler! (l2004h and 
iLenorzer et al.l (120041) . the low value of hv/kT in the IR leads 
to the fact that even small departures from LTE become substan- 
tially amplified (in contrast to the situation in the UV and opti- 
cal). This can be immediately seen from the line source-function 
in the Rayleigh- Jeans limit, 



S L /B v x(l+6/(hv/kT)y 



6 = bilb u - 1 



(2) 



where b\ and b u are the NLTE departure coefficients for the lower 
and upper level, respectively. For temperatures at 30 kK, the 
value of hv/kT is 0.24 at Br y and 0.11 at Biv Thus, under typ- 
ical thin-wind conditions (where in the line forming region the 
lower level becomes underpopulated compared to the upper one, 
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see below) the line source-function can easily exceed the contin- 
uum, or can become even "negative", i.e., dominated by induced 
emission. E.g., for the case of bi/b u = 0.95, Sl/B v « 1.83 at 
Br a , whereas for a ratio of 0.9 already a value of S l/B v « 11 is 
present, and lasering sets in at a ratio of 0.89. 

From these examples, it is immediately clear that the syn- 
thesized profile, particularly Br a , reacts very sensitively to this 
ratio, where the major effect regards the height of the narrow 
emission peak that we will use to constrain the mass-loss rate. 
Thus, we must also check the influence of uncertainties in atomic 
data and atmospheric parameters that can influence this ratio and 
might weaken our conclusions. 

To investigate the general formation mechanism and the 
above problems, we have calculated a large number of mod- 
els exploring the sensitivity of Bi a on various effects, and will 
comment on those in the following, by means of our model of 
HD 37468 with atmospheric parameters as outlined in Table [2] 

General behavior. In Fig. [12] we compare the reaction of U a , 
Br ff , Br y , and Pf y on different mass-loss rates, varied within 
M = 5-10" 10 (solid black) and M0" 7 M o /yr (long-dashed, 
blue). Whereas a clear reaction of H ff is found only for M > 
5-10~ 8 M o /yr, Br„. remains sensitive at even the lowermost val- 
ues. For increasing mass-loss, the height of the emission peak 
in Br ff decreases, whereas the wings (in absorption for lowest 
M) become more and more refilled, going into emission around 
10~ 7 M©/yr. From Fig. [12] we note that also the Br y and Pf y line 
profiles are more sensitive than H a in the thin wind regime as 
their wings typically require a factor of 5 lower M to start dis- 
playing reactions to mass-loss. 

The refilling of the wings with increasing mass-loss can 
be explained by the increasing influence of the bound-free and 
free-free continuum (i.e., the typical continuum excess in stel- 
lar winds becomes visible) as well as a certain "conventional" 
wind emission. However, at line center the line processes al- 
ways dominate and the peak height depends on the location (in 
t) where the wind sets in. This is shown in Fig.[13]which relates 
the strength of the Br^ source function at each point of the atmo- 
sphere with the gradient of velocity field and reveals whether the 
photosphere, transition region or wind control the resulting Br ff 
line profile. 

Depopulation of the n = 4 level in th e outer photosphere. 
As outlined in the introduction, already fAuer &Mihalasl(ll969h 
found in one of their first NLTE-models a stronger depopulation 
of n — 4, compared to n = 5. They argued as follows: Whenever 
the density becomes so low that collisional coupling plays no 
role, the decisive processes are recombination and cascading, as 
in nebulae. Because the decay channel 4 — » 3 is very efficient in 
this process, level 4 becomes stronger depopulated than level 5, 
and the core of Br„ goes into emission. 

Since their findings refer to different conditions (T e g = 
15,000 K), and since it is important to understand the depen- 
dence of the depopulation on the precision of the involved pro- 
cesses, we investigated the process in more detail. As it turned 
out that line-blocking/blanketing effects have a minor influence 
on the principal results for hydrogen (only the absolute peak 
height of Br ff is affected, but not its systematic behavior), we 
used a pure hydrogen/helium model (with parameters as derived 
for HD 37468, but with a very weak wind, M= 5 • 10~ 12 M Q /yr) 
for this purpose, in order to allow for a multitude of calculations. 
The results of our investigation are displayed in Fig. [14] In the 




-2 -3 -4 -5 -6 

Log iau_Ross 



Fig. 13. Line source function of Br„ (in units of the local Planck 
function) and velocity gradient (vertically scaled) as a function 
of tr oss , for the same M-sequence as in Fig. Q~2] Note that the 
transition region between photosphere and wind moves towards 
lower tr oss for decreasing M, i.e., the model with the lowest 
mass-loss is located at the right of the figure. 



optically thin part of the photosphere (10~ 4 < tr oss < 10~ 8 ), 
all departures remain roughly constant, where the ground-state 
(dashed-dotted) is overpopulated by a factor of 10, the 2nd level 
(solid) is roughly at LTE and, indeed, b\ (dashed) is smaller than 
bs (dotted). In contrast, the groundstate is strongly overpopu- 
lated in the wind (oc r 2 for constant temperature, see below), 
whereas the excited levels are overpopulated by factors between 
10 and 5, with a different order than in the photosphere, i.e., 
b 5 < b 4 < b 2 . 

In the following, we will investigate the influence of vari- 
ous effects that determine this pattern, by solving the NLTE rate 
equations and omitting certain rates, but using a fixed radiation 
field (which is legitimate at least in the optically thin part of the 
atmosphere). 

At first, we checked the influence of the collisional rates, by 
leaving them out everywhere. The result of this simulation is 
displayed on the left of Fig. [14] Obviously, below log tr oss = -5 
collisions do not play any role, since the departure coefficients 
with (in black) and without (in color) collisions are identical. In 
the lower photosphere, of course all departures are affected by 
collisions, but particularly for level 4 and 5 a difference is visible 
until log tross = -5, where these levels remain thermalized until 
log tr oss = -2. Since the major part of the line core forming 
region of Br ff is just in the range between 10~ 2 < tr oss < 10~ 4 , 
we have to conclude that collisions do play a certain role in the 
formation of Br ff , and will return to this point later on. 

Before doing so, however, we try to understand the depop- 
ulating processes in the outer photosphere (opposed to the con- 
ditions in the wind) by concentrating on that part which is not 
affected by collisions, in order to gain insight into the dominat- 
ing mechanism. To this end, we display the results of various 
approximations on the right of Fig. [14] by considering the outer- 
most photosphere (indicated by dotted lines on the left panel). 

The first approxim ation ("2") follows the suggestion by 
Auer & Mihalas (1969), i.e., we solved the rate-equations for 
a Case A nebula approximation, i.e., allowed for ground-state 
ionization, spontaneous decays and radiative recombination into 
the excited levels. In this case, we derive departures with 1 > 
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Fig. 14. Underpopulation of the hydrogen n — 4 level (compared to n — 5) in the outer photospheres of late-type O-dwarfs with very 
thin winds (pure H/He model atmosphere). 

Left: Departure coefficients for n = 1,2,4,5 (dashed-dotted, solid, dotted, dashed), accounting for all processes (radiative and 
collisional). Overplotted are the corresponding departure coefficients (red, blue, green, magenta) as resulting from a NLTE solution 
discarding the collisional processes. Note the underpopulation of n = 4 compared to n — 5 in the outer photosphere (responsible for 
the line core emission in Br ff ), which is no longer present in the wind. 

Right: Departure coefficients of n = 2,4,5 (solid, dotted, dashed) in the outermost photosphere (corresponding to the region 
embraced by dotted lines on the left), as resulting from the complete NLTE solution and various approximations, the latter all 
without collisions: 2: nebula approximation (Case A); 3: as 2, but including excitation/induced deexcitation from resonance lines 
(roughly Case B); 4: as 3, but including ionization from excited states; 5: as 4, but including excitation/deexcitation from all lines 
with lower level n < 3; 6: as 4, but including excitation/deexcitation from all lines with lower level n < 5 (see text). 
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bs > b4 > £>2 (i-e., level 4 is more strongly depopulated than 
level 5), but far away from the complete solution. Moreover, the 
ground-state (not displayed) is roughly consistent with the exact 
solution, whereas the departures of the excited levels in the wind 
do not differ from the conditions in the outer photosphere. 

In simulation 3, we switched on the resonance lines, by 
including the corresponding excitation and induced deexcita- 
tion rates (roughly corresponding to Case B nebula conditions). 
Immediately, all departures in the wind obtain values very close 
to the complete model (i.e., bi > bi\ > b$ > 1), whereas in the 
outer photosphere the "correct" order is achieved (Z?2 > b^ > b\), 
though at a much too high level. 

This is cured by simulation 4, where ionization from the ex- 
cited states is switched on. By this process, all levels are depop- 
ulated again, and the corresponding solution looks very close to 
the exact one. (The departures in the wind remain unaffected, 
since these rates are almost unimportant because of the strongly 
diluted radiation field). 

One might conclude now that the remaining missing rates 
(excitation/deexcitation between excited levels) are negligible. 
Unfortunately, this is not the case, at least if one is interested 
in the precise ratio of b^/b^, which is of major importance in 
our investigation. Though the line transitions between the ex- 
cited states are optically thin, the mean line intensity is still large 
enough to be of influence. This becomes clear by comparing 
simulation 5 with 6. In the former, we have included the exci- 
tation/deexcitation from all lines with lower level n < 3 (i.e., 
level 4 and 5 are only affected by line transitions in terms of res- 
onance lines and spontaneous decay), with the effect that now 
b\ ~Z b$, i.e., the emission core of Br ff would disappear. Only if 
all bound-bound processes for lines with at least a lower level of 
n < 5 are included (simulation 6), the "exact result" is recovered 
(i.e., higher levels contribute indeed almost only by spontaneous 
decay). 

In summary, even in the outermost photosphere the actual 
value of Z?5 /&4 is controlled by (almost) all radiative processes 
(with additional collisional contributions in the lower photo- 
sphere), and thus depends on a precise description of the con- 
tinuum and line radiation field and to a lesser extent on the cor- 
rect run of the temperature stratification (entering the recombina- 
tion coefficients). Interestingly, however, we have also seen that 
in almost all of our simulations level 4 is more underpopulated 
than level 5, independent of the various processes considered. In 
Appendix iBl we will show that at least this principal behavior 
can be regarded as the consequence of a typical nebula-like situ- 
ation, namely as due to the competition between recombination 
and downwards transitions. 

In the wind, on the other hand, a Case B like nebula approx- 
imation is able to explain the run of all hydrogen occupation 
numbers alone. As it will be shown in the appendix, the actual 
conditions in the outer wind depend strongly on whether line- 
blocking/-blanketing is considered or not. In all cases, however, 
the abrupt decrease of the Br ff source function in the transition 
region between photosphere and wind is triggered by the onset 
of dilution and the Doppler-effect shifting the (resonance-)lines 
into the neighboring continuum, thus effectively pumping the ex- 
cited levels. 

5.2. Influence of various parameters 

Collision strengths. As already outlined, collisions do play a 
role in the formation of the emission peak of Br^ and, even 
more, in the line wings. In particular, these are the collisional 
ionization/recombination processes for n > 3 and the collisional 
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Fig. 16. Influence of different microturbulent velocities on Br a 
(left) and Pf y (right). See text. 



excitation/deexcitation processes within transitions i <-> j, i > 
3, j > i + 1, (strongest for 7 = 1 + 1) which keep the occupation 
numbers for levels > 3 and higher in of close to LTE than ex- 
pected from considering the radiative processes (for given radia- 
tion field) alone. From comparing the departure coefficients cal- 
culated with and without collisional rates, one might predict that 
a decrease of the collisional strengths in complete models will 
increase the strengths of the absorptions wings (since, around 
tross = 1CT 2 , b\ ~ bs ~ 1 compared to b\ > b$ > 1 with and 
without collisions, respectively), whereas the influence on the 
emission is difficult to estimate. In Fig. [15] we display the result 
of three calculations using different sets of collision strengths 
and its impact on the H a and Br„ profiles. 

Our updated cmfgen "standard" model u tilizes the hydrogen 
collision strengths from Mih alas et al.l (11975 ) (MHA), which are 



compared to our previous data set from lGiovanardi et al.l (Tl 987 



17) 
ler 



(GNP) and recent collision strengths from Przvbilla & Butlei 
(120041^1 (PB). The difference between these data sets is signif- 
icant, particularly for transitions with intermediate i, j such as 
Br a . Our "standard" MHA collision strengths lie in between the 
GNP and PB data sets, the latter being typically a factor of five 
smaller than the GNP implementation. However, the reaction of 
the departure coefficients seems to be small. The "only" differ- 
ence is a weak increase of b\ in the lower part of the line-forming 
region (as predicted above) and a weak decrease of b$ in the 
outer one, i.e., the NLTE effects become increased everywhere. 
Consequently, the absorption wings of Br ff become deeper and 
the emission-peak higher, when using the data set w ith reduced 
collision strengths (PB) bv lPrzvbilla & Butlerl d2004O (Fig. [151 
right panel, dotted profile), where the small differences in depar- 
ture coefficients are (non-linearly) amplified according to Eq. [2] 
Since the changes in collision strengths affect both the peak of 
the line core and the absorption wings of Br a , they cannot be 
mapped directly onto M variations, as the latter only modifies 
the line core in the thin wind case. Such changes may rather be 
accomplished by slightly modifying the gravity of the star. 

Microturbulence. One of the basic "unknowns" in the cal- 
culation of synthetic profiles based o n model atmospheres is 
the microturbulent velocity, Vfnrh (e.g., Smith & Howarth 1998. 
Villamariz & Herrerol 12000. Rep olust et all 120041 iHunter et al.l 
2007). Though it is possible to obtain a "compromise" estimate 
for this quantity from a simultaneous fit to a multitude of differ- 



10 based on ab-initio calculations by Keith Butler 

11 Note that a similar investigation performed by these authors gave 
different results, due to numerical problems (N. Przybilla, priv. comm.). 
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ent lines, it is well known that different lines indicate different 
values, pointing to a dependence on atmospheric height. As a 
rule of thumb, in the parameter range considered here a value of 
Vturb on the order of 10 km s _I seems to be consistent with a vari- 
ety of investigations. Fig.[l6]compares the influence of this quan- 
tity on the synthetic Br ff profile, again by means of our model of 
HD 37468, and a typical mass-loss rate. Though a small impact 
is visible in the blue wing (due to the reaction of the Hei compo- 
nent), the major effect concerns the height of the emission peak, 
which increases for increasing Vturb (varied between 7.5 and 17.5 
km s" 1 ). By comparing with Fig. [12] we see that for (very) thin 
winds (with the line wings well in absorption), an uncertainty in 
Vturb °f ±5 km s _1 (which is a typical value) can easily induce 
uncertainties of a factor of two in the deduced M. We note, how- 
ever, that unlike M which changes "only" the height of the emis- 
sion peak (see Fig.[T2l. microturbulence modifies both its height 
and width (see Fig.fToT-left). A similar effect as for Br„ but with 
lower amplitude is found for Pf y (see Fig. [ToT-right). Therefore, 
provided the spectral resolution is high enough, in principle one 
could separate the effects of M and microturbulence and reduce 
the uncertainty in the final M estimate. 

Most importantly, micro-turbulence affects level n= 4 (sim- 
ilar to the influence of the collisional strengths): The higher 
the turbulence, the sooner and the more effective this level be- 
comes depopulated, increasing the line source-function. Thus, it 
is important to adapt v tul -b already in the atmospheric model (— > 
changes in the occupation numbers), and not only in the formal 
integral, as it is often done with respect to metallic lines and Hei. 

/3-law. Here we concentrate on the effects of the steepness of the 
wind velocity law (expressed in terms of ft) on the line profiles 
from thin wind objects. At first, for very thin (= weak) winds 
such as the one from HD 37468, there is almost no reaction at all, 
since the profile, particularly the line core, is not formed in the 
wind. For winds with a somewhat higher density (as HD 76341) 
where the line core of H [V already reacts, the situation is some- 
what different. Figure[T7]shows corresponding L-Band hydrogen 
lines together with H„, where both ft and M have been modified 
in opposite directions to keep the H„ core at the same depth. 
(For this model, such a combination preserves the Br r core as 
well). Interestingly, however, the cores of Br„ and Pf r still react 
strongly, showing that these lines, together with H a and/or Br r , 
can be confidently used to constrain both ft and M in objects with 
(not too) thin winds and to break the classical M-ft degeneracy. 
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Fig. 17. Influence of ft in thin wind objects (here: our model of 
HD 76341). Three combinations of ft and M which keep the H a 
core at the same depth. The different sensitivity of the cores of 
Br a and Pf r to ft breaks the well-known M-ft degeneracy. See 
text. 



Transition Velocity. A variation of the transition velocity (which 
defines the transition from photosphere to wind) in the thin/weak 
wind case shows only minor effects (not displayed here) in the 
emission core of Bi a , provided such transition takes places at 
reasonable velocities (roughly between 0.05 . . . 0.5 V soun d). For 
the rest of the diagnostic lines considered here, we find no effect 
at all. If the transition is moved deeper into the photosphere, 
i.e., below 0.05 V soun d, the Br^ core starts to display a moderate 
sensitivity. 

Finally, we investigated how Br^ and Pf r respond to clump- 
ing in the thin wind regime. Similarly to the dependence on the 
transition velocity, we found no effects, unless severe clumping 
is already present at the base of the photosphere, in which case 
minor effects start to appear in the core of Br ff . 
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6. Discussion 

6.1. On the reliability of the derived mass-loss rates. 

Since the number of M determinations for OB-stars has signif- 
icantly increased during the last decade, and largely differing 
values even for the same objects can be found in the literature, it 
might be necessary to comment on the reliability of the data pro- 
vided here. Actually, the major origin of differences in the mass- 
loss rates bases on different assumptions regarding the clump- 
ing properties of the wind. These range from unclumped media 
over microclumping (constant or stratified, with different pre- 
scriptions on the clumping-law) to the inclusion of the effects 
from macroclumping and vorosity. 

From the comparison provided in Appendix [A] it becomes 
obvious that in most cases the basic quantity which can be de- 
duced from the observations, namely the optical depth invariant 
Q, is rather consistent within a variety of studies, at least when 
concentrating on p 2 -dependent diagnostics. Diagnostics relying 
mostly on UV resonance lines are more strongly affected, as 
can be seen, e.g., from the differences in the corr esponding Q leti 
value derived here and by iFullerton et al.l d2006l) . who studied 
the Pv resonance line alone. Reasons for such discrepancy are (i) 
the influence of X-ray/EUV emission on the ionization equilib- 
rium, particularly in the mid and outer wind (see below), and (ii) 
the strong impact of macroclumping/vorosity o n the formation 
of resonance lines (Sundqvis t et al.ll2010ll201 fl and references 
therein). 

In this work, we have derived the stellar and wind parame- 
ters from a consistent, multi-wavelength analysis, allowing for a 
rather general clumping law. Since we obtained almost perfect 
simulations of the observed energy distributions, from the UV to 
the IR (and sometimes even the radio regime), we are quite con- 
fident on the quality of the provided values. Corresponding error 
estimates have been quoted in Sect. [3] Even admitting that the 
neglect of macroclumping/vorosity might influence the UV reso- 
nance lines (particularly those of intermediate strength) to a cer- 
tain degree, and that our treatment of the X-ray emission is only 
a first approximation, the ubiquity of excellent fits to features 
from a variety of elements forming in different layers cannot be 
considered as pure coincidence. We are optimistic that problems 
within individual features (which can be disastrous in analyses 
concentrating on such features alone) have, if at all, only a mild 
impact in a multi-wavelength study as performed here. 

One specific problem within our analysis which cannot be 
neglected is the problem with the cores of the optical hydrogen 
and Hel lines, encountered for CygOB2 #7 and #8A as well 
as for a Cam (but see our corresponding comment in Sect. 0). 
Briefly repeated, the problem reflects the fact that within our 
analysis we were not able to obtain a simultaneous fit for both 
the cores of these photospheric features and the IR-lines formed 
in the lower/mid wind. For a perfect fit of the IR features, a low 
M in parallel with a quite large clumping factor (CL; « 0.01) 
is needed, whereas a fit of the photospheric line cores requires 
a larger M accompanied with moderate clumping (CL; a 0.1). 
Further tests are certainly necessary to clarify this problem (i.e., 
the shape of the assumed clumping law might still not be opti- 
mum). In the meanwhile we suggest that the mass-loss rates of 
the problematic objects should be considered as a lower limit, 
and might need to be increased in future work. 

Thin and weak winds. As has been outlined in Sect.Q] the vast 
majority of mass-loss determinations for weak winds and weak- 
wind candidates has been performed by analyzing UV resonance 



lines, since H„ is no longer usable at low M, and the IR has not 
been invoked until now. 

A crucial point regarding the reliability of UV M- 
de terminations for w eak-winded stars has been already provided 
by iPuls et al.l d2008l see their Fig. 20). They calculated the di- 
agnostic UV and H a profiles for a set of thin and weak wind 
models where the mass-loss rate was varied by almost two or- 
ders of magnitude whilst the X-ray luminosity was increased in 
parallel to keep the ionization structure and thus the UV lines 
at the observed level. The changes imposed on the wind did not 
reach the photospheric levels, and thus the UV iron forest did 
not change. Most alarming, however, is their finding that essen- 
tially all profiles could be equally well reproduced with any M 
combined with an appropriate X-ray luminosity, L x . Thus, no 
independent UV mass-loss determination is feasible unless the 
X-ray properties of the star are accurately known. 

This problem, of course, is also present in our analysis. But 
here, we have utilized the Br,,, line as the primary M indica- 
tor, due to its high sensitivity on M, being formed in the up- 
per photosphere for weak winded stars. A similar investigation 
of the same M-L x combinations as for the UV, but using Br,,, 
lead to very promising result (see IPuls et all 12008, their Fig. 
22): In the case of very low M, i.e., a deep-seated line for- 
mation region, Br„ turned out to be basically unaffected by X- 
rays. Even with increasing mass-loss rate, the hydrogen core of 
Br a remains unaffected, at least for canonical X-ray luminosities 
(L x /Lboi < 10~ 4 ). Changes, however, arise for the Hen emis- 
sion component of Br a , caused by the large sensitivity of the 
Hen/Hem ionization equilibrium to X-rays. This problem needs 
to be kept in mind for future analyses. 

Finally, to assess the precision of the derived mass-loss rate 
for our weak-wind object HD 37468, we have to consider the 
major sources of error, as discussed in Sect. 15.21 particularly the 
impact of the hydrogen bound-bound collision strengths. Here 
we estimate a total error of plus/minus 0.5 dex, which is irrele- 
vant at such low mass-loss rates. 

6.2. Stratification of clumping factors 

Since our analysis comprises both p and p 2 diagnostics (with 
all the caveats regarding the potential impact of macroclump- 
ing/vorosity), we are able to provide absolute values for f c \{r) 
as well as for M. T his is qu ite different from (and superior to) 
the investigation by IPUL061 who could derive "only" relative 
values (normalized to an outer wind assumed to be unclumped), 
because of relying on p 2 diagnostics alone. 

Dense winds. For most of our dense wind objects, we found 
rather strong clumping close to the wind-base (which has a sub- 
stantial impact on the actual mass-loss rate, see next subsection), 
whereafter the clumping factor decreases towards the mid/outer 
wind. This finding is similar to the results from lPUL06l at least 
qualitatively. An even quantitative comparison is possible for the 
case of f Pup, the only star among our sample for which we 
could determine the run of the clumping factor througho ut the 
entire wind (Fig. [T8l. When normalizing the results by IPUL061 
to a similar value in the outer wind (/ c j « 5), the agreement with 
our results is excellent and re-assuring, since both investigations 
are completely independent from each other and rely on con- 
siderably different methods. Likewise re-assuring is the fact that 
th e derived stratification o f fair) is very similar to recent results 
bv lSundqvist et al.1 d201 ll) . who performed a consistent analysis 
of Pv and hydrogen/helium recombination lines in the 06I(n)f 
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Fig. 18. Radial stratification of the clumping factor, / e /, for £ 
Pup. Black solid: clumping law der ived from our model f i ts. Red 
solid: Theoretical predictions by iRunacres & OwockH ([2002) 
from hydrodynamical models, with self-excited line dr iven insta- 
bility. Dashed: Average clumping factors derived by iPuls et all 
(120061) assuming an outer wind matching the theoretical pre- 
dictions. Magenta solid: run of the velocity field in units of 
100 km s" 1 . See also Sect. @] 



supergiant A Cep (a cooler counterpart of £ Pup), including the 
consideration of macroclumping and vorosity effects. Also for 
this object, it turned out that clumping peaks close to the wind 
base, with a maximum value of f\ as 28, which is rather close to 

the value derived here for £ Pup (see Fig[T8l. 

Ac c ounting as well for the p revious work bvlCrowfher et alj 
(120021) . iHillier et all (120031) and lBouret et alJ (I2003L l2005h fsee 
Sect. Q]), there is now overwhelming evidence for the presence 
of highly clumped material close to the wind base. These find- 
ings are in stark contrast to theoretical expectations resulting 
from radiation hydrodynamic simulations (Sect. Q]), which pre- 
dict a rather shallow increase of the clumping factor, due to the 
strong damping of the line instabi lity in the lower wind caused 
by the so-called line-drag effect dLucvlll984llOwocki & Rvbickil 
1985). In Fig[18]we compare our present results and those from 
PUL06I with prototypical predictions from radiation hydrody- 
namic simulations, here from the work by IRunacres & Owockil 
(2002). Though there is a fair agreement for the intermediate and 
outer wind, the disagreement in the lo wer wind is striking. A 
similar conclusion has been reached by Sundqvist et all d201 lh . 
and further progress on the hydrodynamic modeling seems to be 
necessary to understand this problem. 



Thin winds. Except for HD 217086, all our sample stars with 
thin winds (i.e., H a in absorption) require no clumping to yield 
excellent fits, and for the former object the derived clumping is 
less (/ v = 0.1) than for the typical dense win d case. A gain, this 
finding is in agreement with the results from PUL06, who con- 
strained the clumping factors for thin winded objects to be sim- 
ilar in the inner and outer wind. Thus and in connection with 
our present results, one might tentatively conclude that most 
thin winds are rather unclumped, which immediately raises the 
question about the difference in the underlying physics. Might 
it be that the wind-instability is much stronger in objects with 
dense winds, e.g., due to (stronger) non-radial pulsations, or is 
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Fig. 19. Wind-momentum rates for our sample stars, as a func- 
tion of log(L/L©). Asterisks correspond to supergiants, crosses 
to other luminosity classes. The two weak wind candidates 
HD 76341 and HD 37468 are denoted by large symbols. The 
solid line s hows the theoretical predictions for the WLR by 
IVink et all ((2000), whilst the dashed one has been shifted to 
lower values by 0.55 dex (see text). 



there a connection w ith sub-surface convection, as suggested by 
ICantielloet al"](l2009h ? 



6.3. Wind-momentum rates 

In Fig. [19] we plot the wind-momentum luminosity relation 
(WLR) for our sample stars, i.e., the wind-momentum rates, 
modified by the factor {R t /R Q ) Q - 5 , as a function of luminosity. 
The results of our analysis suggest a well defined relation, if we 
exclude three outliers, the binary Cyg OB2 #8 A well above and 
HD 7634 1/HD 37468 well below the average relation. Of course, 
a much larger sample needs to be investigated before a final con- 
clusion can be drawn. 

From a linear regression to our results, we obtain an 'ob- 
served' re l ation in parallel to the theoretical predictions by 
IVink et all (l2000L solid), but 0.55 dex (factor of 3.5) lower. Due 
to its large deviation (2 dex!), HD 37468 is certainly a weak 
winded star, whereas HD 76341 (with a deviation of 0.9 dex) 
might be considered as a weak-wind candidate, interestingly a 
lb supergiantQ 

On the assumption of an unclumped outer wind, IPUL06I 
found a good agreement between their "observed" and the the- 
oretical WLR. These results could be unified with ours on the 
hypothesis that the outer regions of most winds are actually 
clumped, with a typical clumping factor on the order of 10. With 
regard to our analysis of £ Pup (and also the hydro simulations), 
this seems to be a reasonable, though somewhat too large value. 

A down-scaling of theoretical wind-momenta as obtained 
he re is also consistent w ith results from the analysis of A Cep 
by ISundqvist et alj d201 lh . who derived a mass-los s rate b eing 
a factor of two lower than predicted by IVink et alJ d2000h . Let 
us finally note that an independent "measurement" of the mass- 
loss ra te of ( Pup based on X-ray line emission by Cohe n et alj 
(2010) resulted in a value of 3.5 TO 6 M /yr, with a lower limit 



12 So far, only dwarfs and few giants have been suggested as weak- 
winded stars. 
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of 2.0 -10 6 M /yr, if the abundance pattern would be solar 
(which is rather unlikely). 

7. Summary, conclusions and future perspectives 

In this pilot study, we have investigated the diagnostic potential 
of L-band spectroscopy to provide strong constraints on hot star 
winds, with particular emphasis on the determination of their 
clumping properties and (actual) mass-loss rates, even for ob- 
jects with very thin (=weak) winds. To this end, we have se- 
cured L ( ' band spectra (featuring Br ff , Pf r and Hei 3.703) for a 
sample of ten 0-/early B-type stars, by means of ISAAC@VLT 
and SpeX@IRTF. The sample has been designed in such a way 
as to cover objects with both dense and thin winds, with the ad- 
ditional requirement that spectroscopic data in the UV, optical 
and H/K band as well as radio observations are present and that 
most of the objects have been previously analyzed by means of 
quantitative spectroscopy. 

For all stars, we performed a consistent multi-wavelength 
NLTE analysis by means of cmfgen, using the complete spectral 
information including our new L-band data. We assumed a mi- 
croclumped wind, with a rather universal clumping law based on 
four parameters to be fitted simultaneously with the other stellar 
and wind parameters. Moreover, we accounted for rotational and 
macroturbulent broadening in parallel. 

For the objects with dense winds, we were able to derive 
absolute values for mass-l oss rates and clumping factors (and 
not only relative ones as in lPULOo) . where Pf r and Bi a proved 
to be invaluable tools to derive the clumping properties in the 
inner and mid wind, respectively. 

For our objects with thin winds, on the other hand, the nar- 
row emission core of Br„(in combination with its line wings) 
proved as a powerful mass-loss indicator, due to its strong re- 
action even at lowest M, its independence of X-ray emission, 
and its only moderate contamination by additional effects such 
as atomic data, microturbulence and velocity law. 

Contrasted to what might be expected, the height of the Br ff 
line core increases with decreasing mass-loss. This is a conse- 
quence of the transition region between photosphere and wind 
"moving" towards lower tr oss for decreasing M, so that more 
and more of the strong line source function becomes "visible" 
when the wind becomes thinner. The origin of such a strong 
source function is a combination of various effects in the up- 
per photosphere and the transonic region that depopulates the 
lower level of Br^, n = 4, stronger than the upper one, n — 5. By 
detailed simulations, we explained this in terms of a nebula-like 
situation, due to the competition between recombinations and 
downwards transitions, where the stronger decay from n — 4 
(compared to the decay from n — 5) is decisive. Pi 

As an interesting by-product, the specific sensitivity of Br (! , 
and Pf r on the velocity field exponent fi in (not too) thin winds 
allows a break in the well-known M-/3 degeneracy when us- 
ing H„ alone. On the other hand, the dependence of Br ff on /3 
for weak winds is negligible, which decreases the error bars of 
the derived M. A major point controlling the depopulation of 
the lower level of Br a and thus the height of the emission peak 
are the hydrogen bound-bou nd collisional data which are used. 
Our models utilize data by Mihala s"et al.l (1 19751) . which provide 
something of a compromise regarding collisional strengths when 
comparing with other data sets. Overall, we estimate the error in 
our Br^-M determination of weak-winded stars by + 0.5 dex. 

13 though most other processes play an important role as well, by es- 
tablishing the degree of depopulation, which controls the peak height. 



This is certainly better than UV diagnostics that strongly depend 
on the assumed description of X-ray properties and are further 
hampered by the impact of macroclumping and vorosity. 

We compared the results of our analysis with those from pre- 
vious work, in particular the derived L e ff, \ogg and Q-values 
(the optical depth invariant(s)). With respect to L e ff, the (aver- 
age) agreement is significantly hampered due to sizable differ- 
ences for Cyg OB2 #8C (AL eff ~ 4000 K), whereas the rest 
agrees within the conventional errors of + 1000 K. The large 
disagreement for the former object has been attributed to differ- 
ences in the derived rotational/macro turbulent velocities. On the 
other hand, the agreement with respect to Q is satisfactory. 

Almost all of our dense wind objects require large clump- 
ing close to the wind base, whereas for the thin winded stars we 
did not need to invoke clumping at all. Our clumping factors for 
the best studied object, £ Pup, agree very well with the work by 
IPUL061 if the clumping in the outer wind (which could not de- 
rived by the latter authors) is scaled to similar values. Moreover, 
our results on strong clumping in the lower wi nd are also con- 
sistent with other findings, particularly those by Sundavist et al.l 
(1201 ll) . and challenge present radiation hydrodynamic simula- 
tions which predict a much shallower increase of the clumping 
factors. 

Because of using Bi a as a mass-loss indicator, we were 
able to fully characterize one weak-winded star, the 09.5 dwarf 
HD 37468, and one weak wind candidate, the 091b supergiant(!) 
HD 76341. 

Finally, our results suggest a well defined WLR (discarding 
the two weak- winded objects and the binary Cyg OB2 #8A) tha t 
is located 0.55 dex below the pred ictions by I Vink et al.l (|2000). 
From a comparison with PUL06, it seems likely that at least 
dense winds are considerably clump ed in the ir outer regions. We 
suggest that the mass-loss rates from lPUL06l are upper limits in- 
deed, and that a downscaling of their values by factors on the 
order of 2 to 3 seems likely. 

From all these results, we conclude that the diagnostic poten- 
tial of IR L-band spectroscopy for deriving clumping properties 
and mass-loss rates of hot star winds is really promising. We 
suggest to extend this rather small sample with further L-band 
observations for a carefully selected sample of OB stars (with 
large and low M) in order to derive statistically conclusive re- 
sults on the "true" mass-loss rates from these stars. Any result 
drawn from only a few objects suffers from a variety of prob- 
lems (e.g., the objects might be peculiar, or problems related to 
the diagnostic tools might remain hidden), and only a careful 
analysis of a large number of objects (performed with the same 
diagnostic tool) allows the identification of trends and outliers. 

Moreover, there is the additional problem of variability. Due 
to its higher sensitivity to mass-loss, it might b e expected that 
Br g i s an even better candidate than H a (e.g., iMarkova et al.l 
120051 and references therein) to study and to analyze wind vari- 
ability in dense winds. The impact of Br a line profile variabil- 
ity on mass-loss/clumping diagnostics needs to be investigated 
as well. To our knowledge, corresponding observations have not 
been carried out so far, and deserve future interest. In thin winds, 
on the other hand, a strong variability is not to be expected, be- 
cause of the near-photospheric origin of Br ff and assuming a sta- 
tionary photosphere. This expectation needs to be confirmed as 
well, and a contradicting outcome might point to variability in 
the location of and the conditions in the transition zone, e.g., 
related to pulsations. 

Finally, and from our experience accumulated so far, we are 
confident that Br„will become the primary diagnostic tool to 
measure very low mass-loss rates at unprecedented accuracy, 
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thus clearly identifying weak-winded stars and quantifying the 
degree with which they lie below theoretical predictions. 
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Appendix A: Detailed comparison of present 
results with other investigations 

The present work comprises a detailed analysis of a small sample of hot stars, 
based on the combination of optical, NIR and UV spectra with one of the most 
sophisticated NLTE atmosphere codes presently available, cmfgen. Thus, a com- 
parison of the derived results with those from more restricted investigations (with 
respect to wavelength range) based on alternative atmosphere codes provides an 
opportunity to address typical uncertainties inherent to the spectroscopic analy- 
sis of such objects, caused by different data-sets and tools. As outlined in Sect.|2] 
(cf. Table most previous investigations of our targets have been analyzed by 
means of fastwindQ or by (quasi-) analytic methods designed for specific di- 
agnostics su ch as Hg or the IR-/mm-/radio continuum (for an overview of these 
methods, see Kudritzki & Puis 2000, Puis et al. 2008 and references therein). 

Brief comments on important differences between our and those results have 
been already given in Sects. |4]and|5] and the complete set of the various stellar 
and wind-parameters is presented in Table lA. 2l In the following, the various in- 
vestigations are referred to following the enumeration provided at the end of this 
table (ref#). Note that a direct comparison of mass-loss rates is still not possi- 
ble, due to the uncertainties in distances (and thus radii) for Galactic objects. 
Moreover, previous analyses were based on either unclumped models (ref# 2-7) 
or that the derived clumping factors had to be normalized to the clumping condi- 
tions in the outermost wind, which are still unclear (ref# 8). Thus, a meaningful 
comparison is possible only for the optical depth invariants F^l 



M 



(R-frVoo) 1 ' 5 V/v 



and 



M 



(A.l) 



which describe the actual measurement quantities related to (i) p 2 -dependent 
processes (Q), if / v is the average clumping factor in the corresponding forma- 
tion region and the clumps are optically thin, and (ii) p-dependent processes 
(Qres), under the assumption that clumping plays only a minor role (but see 
Sect.[TJ. 

Fig. I A. II provides an impression of the differences in the most important 
parameters for the individual stars, by comparing the effective temperatures and 
gravities (upper 3x3 panels), and the optical depth invariants and luminosities 
(lower 3x3 panels). Note that all panels provide identical scales, to enable an 
easy visualization. From the figure, it is quite clear that typical differences in 
T c ff are of the order of 1 ,000 to 2,000 K, with corresponding differences of 0. 1 
to 0.2 dex in log g. The largest differences are found for the effective tempera- 
ture of Cyg OB2 #8C (roughly 4,000 K when comparing with ref# 3 and 6/8), 
which has been already discussed in Sect. [4] and most probably relates to an un- 
derestimation of v sin ; in these investigations. It is reassuring that in most cases 
the connecting lines between our (crosses) and the other results in the T c ff-log g 
plane have a positive slope, indicating that higher temperatures go in line with 
higher gravities and vice versa, which is consistent with the behavior of the grav- 
ity indicators (usually, the wings of the Balmer lines). 

The average differences with respect to effective temperature, (AT c ff), are 
presented in Table lA. 11 discarding ref# 1 and 8 who adopted the stellar parame- 
ters, mostly from ref# 2-7. Major discrepancies seem to be present when compar- 
ing with ref# 3 and 6, who derived temperatures being on average 1 ,700 K higher 
than our results. Note, however, that a large part of this discrepancy is caused by 
the results obtained for Cyg OB2 #8C (see above). The other three investigations 
(ref# 4, 5 and 7) deviate, at least on average, much less from the present one, 
by a few hundreds of Kelvin. The dispersion of the differences, however, is very 
similar in all cases, about 1,700 K, except for ref# 5, with a dispersion of 600 K. 
Thus, overall, the dispersion of (Ar c ff) is somewhat larger than to be expected 
from the typically quoted individual uncertainties of 1 ,000 K, which should give 
rise to a dispersion of 1 ,400 K. 

Regarding the optical depth invariants, the situation is satisfactory. Except 
for ref# 3, the mean differences are at or below 0.1 dex (25%), with a dispersion 
of typicall y 0.3 dex (factor of 2 ), which is consistent with the typical individual 
eiTors (see Markova et al. 2004 for a detailed analysis). This result is particularly 
obvious from Fig. IA.1I (lower panels), where in most cases all investigations 
show rather similar Q values. 

Regarding the Q lcs v alues which are relev ant when comparing with ref# 1 
(the Pv investigation by Fullerton et al. 2006), the discrepancy is still large, 
particularly when accounting for the fact that our results indicate considerable 
clumping, thus reducing the absolute mass-loss rate significantly with respect to 
previous investigations. On average, we find (A log Q m! ) ~-1.4 dex, with a dis- 
persion of 0.8 dex, Thus, either (i) the actual mass-loss rates are even smaller 
than derived here, or (ii) the ionization fraction of Pv (remember that the results 
from Fullerton et al. include the product with this quantity) is very low, of the 



14 which relies on certain approximations mostly related to the treat- 
ment of (EUV-)line-blocking 

15 for a derivation, see Kudritzki '&PuTsll200(]| and lPuls et all2008l 



Table A.l. Mean difference of derived effective temperature, 
Ar e ff = r e ff(ref i) - r e ff(this work) and derived optical depth 
invariant, A log Q, for N objects from reference i,i - 1 , 8 in com- 
mon with our sample (cf. Table lA. 2"l i. tsJ^s in kK, A log Q in dex. 
For ref# 1, we compare the Q KS -values, whereas for the rest we 
compare the g-values normalized by / v (see Eq. IA.ll ) Positive 
values indicate that the results from the specific work are on av- 
erage larger than the results derived here, cr is t he dispersion of 
these differences. No values are given for ref# 2 (Kudritzki et al. 
1999), since there is only one object in common. 



ref# N 


<Ar eff > 


o-(Ar cff ) 


(Alog Q) 


o-(Alogg) 


1 4 






-1.38 


0.76 


3 3 


1.63 


1.72 


-0.21 


0.31 


4 3 


0.13 


1.70 


0.05 


0.28 


5 3 


-0.57 


0.59 


0.11 


0.28 


6 4 


1.75 


1.79 


-0.02 


0.38 


7 8 


-0.21 


1.67 


-0.06* 


0.31* 


8 6 






-0.08 


0.23 



* mean difference and dispersion only for 7 objects, excluding the weak- 
winded object HD 37468. 



order of 4% (which would require extreme conditions in the wind, e.g., a very 
strong X-ray/EUV radiation field), or (iii) the line formation calculations of UV- 
resonance lines (both in the investigation by Fullerton et al. and in our analysis) 
require some additional considerations, such as the presence of opt ically thick 
clumps and/or the inclusion of a porosity in velocity space, see Sect. lll 10 ! 

In summary, we conclude that at least the analysis of Q seems to be well- 
constrained, and that different investigations give rather similar results. The re- 
maining problem is the determination of actual mass-loss rates, which involves 
the "measurement" of (absolute) values of clumping factors. As we have shown 
in this investigation, L-band spectroscopy turns out to be a promising tool for 
this objective. Let us note that only a measurement of actual mass-loss rates will 
enable a strict comparison with theoretical predictions (as performed in Sect.|6j, 
to identify present shortcomings and to provide "hard numbers" for evolutionary 
calculations. 

The precision of effective temperatures, on the other hand, is less satisfac- 
tory. Irrespective of the fact that we did not find a real trend in the average dif- 
ferences with respect to three from five investigations, the dispersion is quite 
large, and individual discrepancies amount to intolerable values. Because of our 
detailed analysis covering a large range of wavelength domains and using a state- 
of-the-art model atmosphere code based on an "exact" treatment of all processes, 
we are quite confident that the 7" c ff -errors in our work are of the order of 1 ,000 K 
or less, which means that the corresponding errors in the previous investigations 
must be of the order of 1,400 K or more. Additionally, two from five investiga- 
tions gave a rather large average difference with respect to our results, which is 
alarming since all five investigations have been performed with the same NLTE 
atmosphere code. Insofar, recent attempts to provide reliable spectral-type-T c ff- 
calibrations have to be augmented by results from large samples to decrease the 
individual scatter in a statistical way. 



Appendix B: Occupation numbers of the hydrogen 
n - 4 and n = 5 level in the outer atmospheres of 
late O-type stars with thin winds 

Conditions in the outer photosphere. As we have seen from Fig. 1141 

almost all of our simulations (and many more which have not been displayed) 
resulted in a stronger depopulation of level 4 compared to level 5 in the outer 
atmosphere, independent of the various processes considered. One might ques- 
tion how far this result can be explained (coincidence or not?). To obtain an 
impression on the relevant physics, we write the rate equations for level i > 1 in 
the following, condensed fornvl again neglecting collisions, and assuming that 
ionization is only possible to the ground-state of the next higher ion (as it is the 
case for hydrogen): 



m 2j AijZij 



Yi n ' A i- 



(B.l) 



16 Another possibility, thou gh less likely, is a strong und erabundance 
of phosphorus, as claimed by Pauldrach et al. ( 1994, 2001). 

17 The case i = 1 will be considered separately below. 
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Fig. A.l. Comparison of the results obtained in this work (crosses) with results from other investigations (for data and reference 
identifiers, see Table IA.2t . Upper 3x3 panels: \ogg vs. T e g; all panels have the same scale, corresponding to an extent of 5000 K 
in T e g and 0.4 dex in logg. Lower 3x3 panels: log Q vs. log(L/L ), with axes extending over 0.4 dex in log(L/L ) and 3.5 dex 
in log Q. In order to facilitate the comparison with p 2 -diagnostics, all Q values have been normalized to fy - 1 (see Eq. lA.U . 
The asterisks provide the gres-values which have to be compared with the corresponding values from ref# 1 dFullerton et al.l f2006. 
p-diagnostics, but including the product with the ionization fraction of Pv). Note that all Q KS values have been scaled by a factor of 
10 to fit into the individual figures. 
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Table A.2. Comparison of stellar and wind parameters as derived in the present analysis with results from previous investigations. 
Units as in Table |2j optical depth invariants, Q and Q res (see Eq. IA.ll ). calculated with M in M yr" 1 , R* in R ^ and in km s _1 . / v 
values from this work as in Tabled i.e., equal to CL1 (Eq.[TJ. f v values from reference (8) (Pui s et al.ll2006l) ) corresponding to the 
clumping factors within the in nermost clumped regi on ("region 2") extending between 1.1 < r/R^ < 2 and assuming an unclumped 
outer wind. For reference (1) dFullerton et a l. 2006)), the quoted results for M and log Q Ki include the product with the ionization 
fraction of Pv. 



star 


ref. 


*- eff 


l°g£ 


R, 




v sin i 


v 




M 


P 


fv 


log(L/L G ) 


log Qres 


log 2 


Cyg OB2 #7 


t.w. 


45.1 


3.75 


14.7 


0.13 


95 


65 


3100 


1.2 


1.05 


0.03 


5.91 


-14.07 


-12.15 




3 


45.5 


3.71 


14.6 


0.30 


105 




3080 


9.86 


0.9 




5.92 




-11.99 




6 


45.8 


3.94 


14.4 


0.21 


105 




3080 


9.98 


0.77 




5.92 




-11.97 




7 


44.0 


3.71 


14.6 


0.10 


145 




3080 


10.0 






5.86 




-11.98 




8 


45.8 


3.94 


15.0 


0.21 


105 




3080 


<4.0 


0.9 


0.2 


5.95 




-12.05 


HD66811 


t.w. 


40.0 


3.63 


18.9 


0.14 


215 


95 


2250. 


2.1 


0.90 


0.03 


5.92 


-13.66 


-11.86 




1 


39.0 




19.4 




219 




2250 


0.44 


0.5 




5.90 


-14.35 






4 


39.2 


3.65 


17.5 


0.20 


203 




2300 


6.4 


0.92 




5.82 




-12.10 




5 


39.0 


3.59 


19.4 


0.20 


220 




2250 


8.8 


0.9 




5.90 




-12.02 




7 


39.0 


3.59 


19.4 


0.17 


220 




2250 


8.77 






5.90 




-12.02 




8 


39.0 


3.6 


18.6 


0.20 


220 




2250 


4.2 


0.70 


0.2 


5.86 




-11.96 


Cyg OB2 #8C 


t.w. 


37.4 


3.61 


14.3 


0.10 


175 


90 


2800 


2.0 


1.30 


0.10 


5.56 


-13.75 


-12.10 




3 


41.0 


3.81 


13.3 


0.09 


145 




2650 


2.25 


0.9 




5.65 




-12.47 




6 


41.8 


3.74 


13.3 


0.13 


145 




2650 


3.37 


0.85 
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7 


39.0 


3.62 


13.3 


0.10 


145 




2650 


2.0 
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-12.52 




8 


41.8 
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0.13 


145 




2650 


<3.5 


1.0 
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37.6 


3.52 
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0.10 


110 
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2700 
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3 


38.5 


3.51 


27.9 


0.10 
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2650 
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6 
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29.0 
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1550 
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1550 
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2 


28.5 


3.00 
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0.1 


80 




1600 


2.40 


1.25 




5.86 
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<29.0 


3.0 


35.0 


0.1 


80 




1600 


5.25 






5.89 




-12.40 


HD 217086 


t.w. 


36.8 


3.83 


8.56 


0.1 
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80 


2510 


0.028 


1.2 


0.10 


5.08 


-15.28 


-13.55 




1 


36.0 




8.6 




332 




2550 


<0.00174 


1.0 




5.05 


-16.51 






5 


36.0 


3.72 


8.6 


0.15 


350 




2550 


<0.23 


0.8 




5.05 




-13.15 




6 


38.1 


4.01 


8.3 


0.09 


350 




2550 


0.21 


1.27 




5.12 




-13.17 




7 


36.0 


3.78 


8.6 


0.15 


350 




2550 


<0.09 






5.05 




-13.56 


HD 36861 


t.w. 


34.5 


3.70 


13.5 


0.11 


45 


80 


2175 


0.28 


1.3 


1.0 


5.37 


-14.36 


-13.25 




1 


33.6 




15.1 




74 




2400 


0.0013 


0.7 




5.42 


-16.83 






4 


33.6 


3.56 


17.2 


0.1 


66 




2400 


0.97 


0.8 




5.53 




-12.94 




8 


33.6 


3.56 


14.4 


0.10 


66 




2400 


<0.4 


0.9 


0.5 


5.38 




-13.06 


HD 76341 


t.w. 
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0.1 


63 


80 


1520 


0.065 


1.2 


1.0 


5.64 


-14.88 


-13.95 


HD 37468 


t.w. 


32.6 


4.19 


7.1 


0.1 


35 


100 
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0.0002 


0.8 


1.0 


4.71 


-16.90 


-15.74 




T 


30.0 


4.0 


7.1 


0.1 


80 




2300 


0.0165 


1.0 




4.57 




-14.10 



References: m lFullerton et alJd2006h : (2) lKudritzki et alj (fl999. unblanketed analysis); (3)lHerrero et alJd2002h: (4)lMarkova et a"fld2004L based 
on ste llar parameters calibrated to th e results from op tical NLTE analyses by Repolust et al. 2004) (5) Repolust et al .1 d2004l) 7(6) Mokiem et al. 
d2005h : (7) lReporust et alj d2005h : (8) lPuls et al] d2006» . 

n) using the high luminosity solution; h) stellar radius and corresponding quantities scaled to the solution by (5,7) to facilitate the comparison; c) 
stellar radius from this work. 



where A;; are the Einstein-coefficients for spontaneous decay, R/i, and the 
rate-coefficients for ionization/recombination, and Zjj the net radiative brackets 

for the considered line transition, Z,, = (1 — (the fraction denotes the ratio of 

J O jj 

mean line intensity and line source function), n* = nkn e <bjk(T e ) denotes the LTE 
population of level (', accoun ting for the actu al electron and ion density, n e and nt 
(for further details see, e.g., Mihalas 1978), such that the departure coefficients 
are given by i>, = nj/nt. Note that for purely spontaneous decays Zjj = 1, for 
lines which are in detailed balance, Zy = 0, and for levels which are strongly 
pumped (e.g., by resonance lines with a significantly overpopulated lower level), 
Zjj < 0. Solving for the departure coefficients, Eq. lB.ll results in 



ZiKjnj/n'AjiZji+Ria 
2 i> i ™i jZjj + Rjk 



(B.2) 



The sum in the nominator corresponds to the net-contribution of lines from 
"above" (i.e., with upper levels j > i), normalized to the LTE population of the 
considered level, whereas the sum in the denominator is the net-contribution of 
lines to lower levels (J < i). The complete fraction can be interpreted as the ratio 
of populating and depopulating rates, which can be split into the contributions 
from bound-bound and bound-free processes, 



£,•<;(...) 



£;>;(... ) + «;* £;>;(...) + «* 



(B.3) 



For all our simulations 1-6 we have now calculated those two terms which de- 
termine &4 and i>5 . At first, let us concentrate on the outer photosphere, as on 
the right of Fig. 1 141 In almost all cases (except for simulation 5), the 2nd term 
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dominates the departure coefficient, and, moreover, the first term (the ratio!, not 
the individual components) remains rather similar, of order 0.2. Consequently, 
the stronger depopulation of level 4 compared to level 5 is due to the fact that the 
quantity 



Rki 



(B.4) 



is usually larger for level 5 than for level 4, even though Ri& < Rm- the accumu- 
lated transition probability from level 4 to lower levels (A41Z41 +A42Z42+A43Z43) 
is much larger than the corresponding quantity from level 5 to lower levels 
(A51Z51 + A52Z52 + ...). This behavior, finally, can be traced down to the run 
of the oscillator-strengths in hydrogen: On the one side, e.g., A41 is larger than 
A51, etc., whereas, on the other, the corresponding net radiative brackets (Z4; vs. 
Z$j etc.) do not differ too much. 

Two examples shall illustrate our findings. For the complete solution, the 
first term in Eq. IB.3l is roughly 0.22, whereas the 2nd term amounts to 0.62 for 
level 5 and to 0.53 for level 4. Thus, b 4 a 0.75 and b$ as 0.84. For simulation 2, 
with Z,j = Zjj = 1 and /?# = for i > 1, the first term x 0. 18, and the 2nd one is 
0.41 and 0.3, respectively, such that b 4 a 0.48 and b s a 0.6 (cf. Fig.QA}. 

In conclusion, the stronger depopulation of level 4 compared to level 5 in 
the outer photospheres of hot stars can indeed be regarded as the consequence 
of a typical nebula-like situation, namely as due to the competition between re- 
combination and downwards transitions. Different approximations regarding the 
contributing lines do control the absolute size of the departures, but not the gen- 
eral trend. 



Conditions in the Wind. Though the formation of the emission peak of 
Br ff for objects with thin winds is controlled by the processes in the upper part 
of the photosphere, it is also important to understand the conditions in the wind, 
since, as we have seen in Fig. 1131 the onset of the wind prohibits a further growth 
of the corresponding source function: Immediately after the transition point be- 
tween photosphere and wind, the source function drops to values corresponding 
to the local Planck-function (i.e., the departure coefficients of n 4 and ns become 
similar). Only in the outer wind the source function increases again (in contrast 
to the predictions of the pure H/He model, see below), which remains invisible in 
the profile, due to very low line optical depths, //'this abrupt decrease would not 
happen, the monotonic behavior of the strength of the emission peak (Fig. 1121 
would no longer be warranted for mass-loss rates at the upper end of the scale 
considered here, and an important aspect of its diagnostic potential would be 
lost. 

Let us first concentrate on the conditions in the outer wind, where the 
ground-state has a major impact. We stress again that we are dealing here with 
(very) weak winds, i.e., the continuum-edges are formed deep in the photo- 
sphere, whereas the wind and the transition region are already optically thin. 
Otherwise, we could no longer assume a "given" radiation temperature in the 
continuum, but would have to account for a simultaneous solution of radiation 
field and occupation numbers, as it was done, e.g., to explain the ground-state 
depopulation of Hen in dense hot star winds by lGabler et alj 1 1989). 

Within our assumptions of ionizations to the ground-state of the next higher 
ion only and neglecting collisional processes, we obtain an alternative formula- 
tion of the rate equation for the ground-state, by summing up the rate equations 
for all levels ;' (Eq. lB. U . 



For non-blocked models (see Fig. 1 141 . the wind-temperature is not too dif- 
ferent from r ra( ], and b\ becomes strongly overpopulated oc l/W = r 2 . Moreover, 
all net radiative brackets coupled to the ground-state, 
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become strongly negative, since (i) b\ is severely overpopulated and (ii) the 
Doppler effect in the wind allows for an illumination by the continuum blue- 
wards from the resonance-line rest-frame frequencies Vn, i.e., / = WBr(7/ m d,jl) 
(optically thin (Sobolev-)approximation), with 7" rad j| 2> 7/ ra[ ] due to missing 
line-blocking. 

The consequence for the population of the excited levels (Eq. IB. 21 is 
twofold. Because of the strong pumping by the resonance lines, the (normalized) 
population of the higher levels (rij/n*,j > i) is much larger than in the photo- 
sphere, and the line term becomes larger than the recombination term. Second, 
the denominator decreases significantly, due to the direct effect of Zj\ and since 
the the ionization rates oc W become negligible. 

In total, now the first term dominates in Eq. IB. 31 and the situation is just 
opposite to the conditions in the outer photosphere: The lower the considered 
level ;', the larger is the nominator and the smaller the denominator, such that we 
obtain the sequence 62 > 63 > £74. . . > 1 (cf. Fig. 1 141 ) 

For line-blocked models, on the other hand, the cooling by the enormous 
number of metallic lines leads to a strong decrease of the electron temperature 
in the outer wind, and T e becomes much smaller than the radiation temperature 
in the Lyman continuum (for our late O-type model, 10,000 K vs. 25,000 K). In 
this case, ionization, though diluted, outweighs recombination (the exponential 
term in Eq. lB. 71 . and the ground-state even becomes underpopulated (b\ — »0.5). 
Consequently, the resonance lines can no longer pump the excited levels (even 
more, since for blocked models the radiation temperatures close to the resonance 
lines, T m d ,j 1 , are much smaller than in the unblocked case). Thus, we find a situa- 
tion similar to that in the outer photosphere, namely that the 2nd term in Eq. lB. 31 
is the decisive one, and n$ > which is obvious also from the final increase of 
the line-source function for Br ff in Fig. ll2l for all mass-loss rates considered. 

Finally, in the region between the outer photosphere and the outer wind, the 
dilution of the radiation field is faster or similar to the decrease of T e , both for 
the blocked and the unblocked models. Thus, the departure coefficients of level 
4 and 5 increase in this region (due to an overpopulated ground-state, and effec- 
tive pumping due to the onset of the Doppler-shift), though at a rather similar 
rate, with 64 S b$ . Consequently, the source-function approaches the LTE level, 
which explains its abrupt decrease in the transition region. 



niRi k + ^ n,R ik = n\R k \ + ^ n*R ki , 
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since the line contributions cancel out. Solving for the ground-state departure 
coefficient, we find 



Zi>i n*(R ki - biR ik )ln\ + R kl 
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which can be approximated by the well known expression 
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where 7" ra( j is the radiation temperature in the ground-state (= Lyman) contin- 
uum, v > vq, and W the dilution factor. A correction factor of order unity ac- 
counts for the ioniz ation/recombination from the excited levels (for details, see, 
e.g. JPuls et al.120051) . 

From here on, we have to divide between line-blocked and unblocked (e.g., 
pure H/He-) models, as they behave different in the wind (though similar in the 
outer photosphere), due to a considerably different run of electron temperature 
and radiation-field on both sides of the Lyman edge. 



